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INTRODUCTION 

Nearly  250  years  ago  Reaumur  (1735)  recognized  that  the  development 
of  most  plants,  animals  and  insects  was  affected  by  various  environmental 
factors  such  as  soil  characteristics,  nutrient  levels  and  atmospheric 
variables.  There  is  an  optimum  set  of  these  environmental  factors 
specific  for  each  living  organism  which  will  produce  the  maximum  rate  of 
development.  The  important  environmental  factors  for  plants  include  such 
elements  of  climate  as  temperature,  humidity,  rainfall,  wind  and  sunshine. 
Other  factors  such  as  soil  nutrient  levels,  diseases,  insects,  interspecies 
competition,  day  length,  etc.  futher  complicate  the  situation. 

If  we  consider  each  specific 
factor  separately,  by  assuming  that  all 
of  the  factors  except  the  one  under 
consideration  are  near  optimum  for 
the  organism  and  thus  not  limiting  to 
growth  and/or  development  and/or 
activity,  we  can  define  certain 
critical  values  of  the  variable  which 
can  be  related  to  the  development 
and/or  activity  (Fig.  1).   In  the 
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Figure  1 
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figure,  A  is  the  lower  limit  of  the  variable  below  which  no  appreciable 
growth  or  activity  is  in  evidence,  B  is  the  optimum  value  at  which  the 
maximum  rate  of  development  or  activity  occurs  and  C  is  the  upper 
limit  beyond  which  death  or  damage  to  the  organism  occurs.  The  shape  of 
the  development  curve  will  vary  depending  on  the  specific  variable  under 
consideration  and  the  organism  in  question. 


TEMPERATURE  RELATIONSHIPS/PLANT/INSECT 

During  the  last  hundred  years  beginning  with  Fritsch  (1861)  and 
Abbe  (1905)  considerable  effort  has  been  made  by  researchers  in  relating 
environmental  temperatures  to  plant  development.  J.  Y.  Wang  (1963) 
reviewed  some  of  the  many  approaches  that  have  been  reported  in  the 
literature.  One  of  the  more  frequently  used  ideas,  probably  due  to  its 
simplicity,  assumes  that,  within  a  range  of  environmental  temperatures 
commonly  found  in  nature,  plant  growth  and/or  development  can  be 
related  to  temperature  quite  satisfactorily  by  a  straight  line,  for 
example  an  extension  of  the  line  AD  in  Fig.  2.  The  accumulation  of 
energy  represented  by  this  line  has  been  given  several  names:  heat 
units  (HU),  energy  units  (EU)  or 
growing  degree  days  (GDD)  being 
among  the  most  common.  This 
relationship  assumes  that 
regardless  of  how  high  the 
temperature  may  rise  the  rate  of 
development  increases  linearly  with 
increasing  temperature.  The 
accumulation  of  GDD  units  for  any 
given  24  hour  period  is  obtained 
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by  subtracting  a  base  temperature  represented  by  TA  for  the  specific  plant 
species  from  the  Average  Temperature  for  tfTe  24  hour  period. 

More  recently  the  growth  curve  has  been  approximated  by  two  straight 
lines  AD  and  DE  (fig.  2)  where  DE  represents  an  average  of  the  higher 
temperatures  commonly  observed.  This  method  of  accumulating  energy  units 
assumes  that  at  temperatures  above  the  values  TD  the  rate  of  development 
remains  constant.  These  GDD  energy  units  are  calculated  by  addition  of 
the  maximum  temperature  for  the  day  and  the  minimum  temperature  and 
dividing  the  sum  by  two  to  obtain  an  average  but  if  the  maximum  temperature 
exceeds  TD  then  the  maximum  temperature  is  set  equal  to  TD  and  if  the 
minimum  is  less  than  TA  it  is  set  equal  to  TA.   The  average  thus  calculated 
is  then  used  in  the  same  manner  as  described  in  the  first  method  of 
sumation.  Comparisons  of  the  two  methods  at  Davis,  California  indicate 
that  the  predictive  accuracy  for  the  date  of  the  first  bloom  of  tomatoes  is 
improved  by  about  a  half  a  day  using  the  second  method  of  calculation. 
Using  the  same  constants  in  differing  climates,  the  improvement  is 
considerably  greater.  (Richardson  et  al .  1980) 

In  their  1980  study,  Richardson  et  al .  compared  the  accuracy  of 
several  methods  of  relating  the  temperature  environment  to  the  growth  and 
development  of  tomatoes.  A  data  set  for  tomatoes  developed  by  Dr.  Floyd 
Ashton  of  Davis,  California  which  consisted  of  seven  plantings  with  9 
plants  in  each  planting  was  used  to  test  the  models. 

The  models  tested  were  the  two  previously  discussed  plus  two 
curvelinear  models  which  represent  versions  of  the  curve  in  figure  1. 
The  authors  proposed  a  cosine  equation  in  the  form  of  equation  1. 

GDH   .J£     {1    +    COS    [tt+it(™    :    %)]}  (1) 
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Where: 

GDH  =  the  accumulation  of  growing  degree  hours  during  an  hour  when: 

TH  =  hourly  temperature 

TB  =  the  base  temperature  below  which  no  appreciable  growth  will 
occur.  This  corresponds  to  point  A  on  the  curve  in  Figure  1 

TU  =  the  optimum  temperature  or  the  temperature  at  which  the  maximum 
rate  of  growth  will  occur  and  corresponds  to  point  B  on  the 
curve  in  Figure  1 

TC  =  the  critical  temperature  or  the  temperature  beyond  which  little 
or  no  growth  will  occur  and  above  which  irreversible  damage 
may  be  done  to  the  plant  or  organism. 

A  =  TU  -  TB  (the  amplitude  of  the  growth  curve) 

F  =  A  stress  factor  which  can  be  used  to  represent  any  form  of. 
plant  stress  such  as  atmospheric  stress  from  low  humidity, 
soil  moisture  stress  j  stress  from  disease,  insect  damage  or 
lack  of  nutrients. 

A  review  of  the  response  of  several  plant  species  to  their  environmental 
temperatures  indicated  that  rather  than  following  the  true  cosine  curve  to 
TC,  the  rate  of  response  to  increasing  temperature  decreased  more 
rapidly  above  the  optimum  temperature  than  the  rate  of  increase  with 
increasing  temperature  below  the  optimum.  Leavitt  (1980).  Also 
hypothesises  a  similar  response  in  his  curve  on  page  361. 

A  second  response  equation  (2)  was  developed  to  account  for  this 
change  in  the  response  curve. 

GDH  =  FA  {  1  +  cos  [  |  +  |  (  ™  I  jjj  )1)     •  "(2) 

If  the  value  of  TH  is  less  than  TU,  then  equation  1  is  used  in  the 
accumulation  of  GDH;  if  the  value  of  TH  is  greater  than  TU,  then  equation 
2  is  used. 

Testing  the  four  methods  of  calculating  GDD  or  GDH  on  the  tomato 
data  set  showed  a  marked  improvement  in  the  predictability  when  using  the 
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combination  of  equations  1  and  2.  The  average  difference  between  the 

predicted  and  observed  dates  of  first  bloom 'ranged  from  4.27  days 

using  the  original  linear  response  to  only  1.81  using  the  combination. 

TABLE  1 

Comparison  of  predicted  and  observed  dates  of  first  bloom 
of  seven  plantings  of  tomatoes  at  Davis,  California 
as  calculated  from  four  different  models 

Cardinal  Average  difference  between 

Temperatures  Predicted  and  Observed  Dates 

of  First  Bloom  of  Tomatoes 


TB      TU      TC 


10°C 





10°C 

30°C 



10°C 

30°C 

50°C 

10°C 

26°C 

36°C 

4.27  days 
4.02  days 
3.65  days 
1.81  days 


Where:  TB  =  the  base  temperature 

TU  =  the  optimum  temperature 
TC  =  the  critical  temperature 
all  as  defined  in  the  text. 

Similar  increases  in  the  accuracy  of  prediction  of  phenological 

stage  of  several  native  weed  species  competitive  with  tomato  were 

found.  It,  therefore,  appears  reasonable  that  if  the  cardinal  temperatures 

for  the  desired  range  species  can  be  determined;  the  use  of  the  model 

which  combines  equation  1  and  2  will  produce  considerable  increase  in 

accuracy. 

PHENOCLIMATOGRAPHY  MODEL  FOR  SUMMER  RANGES  IN  UTAH 


Under  a  grant  from  the  U.S.  Forest  Service,  Dr.  W.  A.  Dugas,  Jr.  (1980) 
developed  a  model  to  predict  range  plant  development  and  production  using 
Equation  1.  This  model  was  developed  and  tested  using  data  from  the 
Ephraim  and  Straight  Canyon  Watersheds  in  Central  Utah.  These  ranges  are 
all  high  elevation  summer  ranges  and  the  percent  of  plant  cover  is  usually 
much  higher  than  is  found  on  the  lower  elevation  BLM  winter  ranges  in  the 
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central  portion  of  the  state. 

The  objective  of  this  study  was  to  develop  mathematical  equations  to 
predict  range  development  and  production  utilizing  only  the  weather 
variables  commonly  measured  at  National  Weather  Service  climate  stations. 
These  variables  are:  daily  maximum  and  minimum  temperature,  precipitation 
and  snowfall.  The  hourly  temperatures  required  by  equation  1  were 
generated  from  a  linear  interpolation  of  successive  daily  maximum  and 
minimum  values  (Richardson  et  al ,  1974  and  Heuer  et  al ,  1978).  This 
method  assumes  that  the  daily  maximum  and  minimum  temperatures  occur  12 
hours  apart  each  day.  The  hourly  temperatures  for  the  first  12  hours  of 
any  day  are  calculated  from  equation  3. 

TH-  TNi  -  fri  -  1)(TX.  -  TNt)  /  11  (3) 

for  values  of  n  =  1  to  n  =  12.  i  =  the  day  number  of  the  period  for  which 
calculations  are  being  made.  Values  of  TH  for  the  second  12  hours  of 
the  day  are  calculated  from  equation  4. 

THn  =  TXi  -  (n-l)(TXi  -  TN.+1 )  /  11  (4) 

again  for  values  of  n  =  1  to  n  =  12. 

Empirical  relations  were  used  to  determine  the  soil  moisture 
compensation  factor  (F  )  Dugas  (1980).  The  daily  soil  moisture  budgetary 
procedure 

ST,  -  ?.   -   ET,  +  ST,  .  ,  (5) 

where  ST  is  the  storage  of  water  within  the  soil  profile  (equivalent  depth), 
P  the  precipitation  and  ET  the  actual  evapotranspiration  for  day  (i). 
Measurements  or  estimates  of  initial  ST,  daily  values  of  P  and  ET,  and 
the  wilting  point  and  field  capacity  values  for  the  soil  at  each  site  are 
required. 

Daily  potential  evapotranspiration  (ETP  was  calculated  from  the 
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method  of  Blaney  and  Criddle  (1962)  and  actual  evapotranspi ration  (ET) 
from  the  method  of  Hanks  (1974).  Initially;  ETP  was  separated  into 
potential  evaporation  (EP)  and  Transpiration  (TP)  by  an  empirical 
relationship  dependent  upon  the  crop  growth  throughout  the  season.  The 
specific  equations  used  in  the  summer  range  program  are  as  follows: 

Calculation  of  Potential  Evapotranspiration  (modified  Blaney-Criddle)(ETP. ) 

TX.  +  TN.  , 
ETP.  -  (  \       1  }   PDm  (6) 

100  [ND  ] 
Calculation  of  Potential  Evaporation  (EP.): 

^i  =  C  cox  +  co2  (Mf )  +  co3  (I5|fi)2  +  ,(7) 


ETP 

<*  (fig*)**  co5  (™|^l  [  ^i  i 

Where:  PDm   =  Percent  of  daylight  hours 

ND   =  Number  of  days  in  the  month 

CO     for  n  =  1  to  5  -   Constants  for  estimating 
potential  evaporation  bond  on  Hank's  model 

Today  =  Todays  day  number  of  growing  season 

NDGS  =  Number  of  days  in  growing  season  at  site 


Calculation  of  Potential  Transpiration  (TP.): 

TPi  =  ETPi  -  EPi     ,  (8) 

Calculation  of  Actual  Evaporation  (EA. ): 

EAi  *  EPi  <re'1/2  '  (9) 


Calculation  of  Actual  Transpiration  (TA. ): 

™i  ■  TPi  •  AW  •  2  "  •        '(10) 


if  (TA,  >  TPJ  TAj  =  TP. 
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Calculation  of  Actual  Evapo trans pi  rati on  (ETA, ): 

ETAi  =  TA.  =  EA. 

Where:  N  Days  =  number  of  days  since  last  precipitation 
totaling  .20  inches  or  more. 


(11) 


Calculation  of  Present  Soil  Moisture  (PSM): 

PSM  =  BSM  -  ETA  +  SP 
if  PSM  >  FC  PSM  =  FC 
if  PSM  <  WP  PSM  =  WP 

Where  BSM  =  Soil  moisture  of  previous  day 
SP  =  Daily  precipitation 

Calculation  of  Percent  of  Water  Available  (AW): 
AW  =  1  =  (FC  -  PSM)  /  (FC  -  WP) 

Calculation  of  the  Value  of  Compensation  Factor  (Fs) 

Fs  =  y  +  [(1  -  y)  /  cw  1  AWi 

Where  y  =  Potential  minimum  value  of  Fs 

and  cw  =  Critical  percentage  of  soil  moisture  below 

which  the  plant  is  subject  to  moisture  stress 

For  the  summer  ranges  Dugas  estimates  the  values  of: 

y  =  0.5  and 
cw  =  .25 


(12) 


(13) 


(14) 


WINTER  RANGE  MODELS 

J 

Weather  conditions  on  the  winter  ranges  of  the  BLM  are  quite 
different  from  those  found  on  the  higher  elevation  summer  ranges. 
Further,  the  plants  commonly  found  on  these  two  ranges  differ  appreciably 
in  their  production  capabilities  and  their  ability  to  withstand  such 
stresses  of  weather  extremes  as  drought  and  high  temperatures. 

Again  the  major  objective  of  this  study  was  to  determine  the 
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feasibility  of  developing  mathematical  models  to  predict  range  development 
and  production  utilizing  only  the  more  common  meteorological  elements 
measured  at  National  Weather  Service  Climate  stations. 

The  studies  on  tomatoes  and  competing  weeds  briefly  reviewed 
previously  indicated  a  potential  for  marked  improvement  in  the  predictive 
capabilities  of  the  models  if  a  combination  of  the  two  equations,  1  and 
2,  were  used  in  place  of  the  single  equation  used  in  the  summer  range 
model . 

To  properly  utilize  these  equations,  information  on  the  three,  cardinal 
temperatures  for  each  species  is  required.  Since  such  information  was  not 
available,  a  thorough  search  of  the  relevant  literature  was  first  made  to 
obtain  the  best  estimate  of  these  temperatures  and  information  on  other 
factors  which  influence  the  growth  and  development  of  the  desired  species.. 
A  brief  review  of  the  results  of  this  literature  search  follows: 
LITERATURE  REVIEW  ON  THE  SELECTED  SPECIES 
Atriplex  'canescens— Four  Wing  Salt  Bush  (Springfield,  1963) 

Atriplex  canescens  is  one  of  the  most  palatable  of  our  southwestern 
shrubs.  The  leaves,  stems  and  fruits  are  cropped  by  almost  all  classes  of 
livestock.  The  plant  has  a  deep  root  system.  Mature  plants  will  reach 
depths  of  5  to  15  meters  in  alluvial  soils.  Its  extensive  root  system 
makes  the  plant  remarkably  drought  resistant.  (Van  Dersal  1938) 
Afterripening 

Seeds  of  four  wing  salt  bush  apparently  undergo  afterripening.  This 
process  involves  chemical  changes  that  remove  germination  blocks.  Seeds 
of  many  species  are  dormant  at  the  time  of  maturity.  Four  wing  salt  bush 
seeds  frequently  germinate  better  the  second  or  third  year  after 
collection.  From  the  description  of  the  afterripening  process,  it  is 
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possible  a  chilling  requirement  may  be  necessary;  but  Springfield  (1963) 

does  not  seem  to  think  so. 

Germination 

The  uptake  of  water,  imbibition,  is  the  first  process  that  occurs 
during  germination.  Studies  indicate  that  15  to  24  hours  soaking  in 
water  is  required  to  complete  this  process  after  which,  if  environmental 
factors  are  correct,  germination  continues. 
Temperature  effects 

Four  wing  salt  bush  usually  germinates  best  at  relatively  low 
temperatures.  The  woody  plant  seed  manual  (Forest  Service  1948)  recommends 
77°F  day  and  50°F  night  temperatures  for  20  to  30  days  to  test 
germination  of  the  seed.  Hervy  (1955)  obtained  47%  germination  of  four 
wing  salt  bush  at  68°F,  45%  at  59°F  and  77%  at  39°F. 

Hilaria  jamesii— Galleta  Grass  (West  et  al  1972) 

Galleta  grass  is  a  normally  course  perennial  grass  varying  in  height 
from  8  to  65  cm  (3  to  20  inches).  One  of  the  most  characteristic  features 
of  the  species  is  its  tough,  woody  rhyzomes  which  represent  the  plant's 
surest  means  of  reproduction.  Most  roots  and  rhyzomes  are  found  above 
50  cm  (18  inches)  in  the  soil.  Roots  commonly  extend  up  to  40  cm  beyond 
the  crown  in  all  directions. 
Germination 

Galleta  seed  collected  from  various  areas  in  the  west  germinated  85% 
or  higher  at  temperatures  ranging  from  16  to  32°C  (60-90°F). 
Vegetative  Growth 

A  study  by  Jameson  (1965)  showed  that  galleta  maintained  a  green  leaf 
height  greater  than  about  5  cm  (2  inches)  from  April  20  to  November  1  in 
northern  Arizona.  The  plants  reached  a  summer  growth  plateau  at  a  height 
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of  12  to  13  cm  (4.5-5.5  inches)  and  reached  an  average  peak  height  of 
38  cm  (15  inches)  about  September  11. 

Jameson  also  reported  that  the  early  summer  decline  in  growth  rate, 
beginning  of  the  summer  growth  plateau  varied  with  elevation  and  probably 
with  temperature.  He  concluded  that  the  time  of  growth  corresponded 
closely  to  the  seasonal  precipitation  pattern,  but  earlier  than  normal 
precipitation  did  not  result  in  earlier  than  normal  growth  of  a 
warm  season  grass  like  galleta.  (Warm  season  disagrees  with  later 
information  reported  on  cardinal  temperatures  in  the  same  bulletin  which 
are:  TB  =  4°C,  TU  =  26°C  and  TC  =  40°C. 

The  grass  undergoes  fall  regrowth  with  adequate  rainfall  after  the 
summer  dormancy  which  appears  to  be  drought  induced.  When  adequate 
moisture  is  available  during  the  late  summer  growing  season,  growth 
continues  throughout  the  season. 

The  time  of  seed  formation  varies  considerably.  In  dryer  areas, 
the  plants  will  produce  seed  and  become  dormant  a  full  month  before  they  do 
in  more  normal  environments.  "Growth  initiation  and  vegetative  production 
appears  to  be  temperature  related,  but  flowering  appears  to  be  day  length 
triggered.  Growth  initiation  appears  to  be  associated  with  average 
maximum  temperatures  in  the  air  of  16°C  (60°F).  The  species  can  withstand 
high  temperatures  and  a  desicating  environment.  The  first  seed  heads 
formed  in  Tucson  on  March  24,  1970,  but  in  the  greenhouse  and  fields 
generally  seed  formation  ranges  between  June  and  October.  Thus  the 
species  evidently  has  some  latitude  in  terms  of  daylength  to  trigger  seed 
formation.  It  appears  to  flower  whenever  and  as  long  as  adequate  moisture 
is  available.  The  lethal  temperatures  are  near  61 °C  (142°F). 
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Oryzopsis  hymenoides — Indian  Rice  Grass 

Oryzopsis  growth  and  production  appears  to  be  determined  more  by 
spring  moisture  and  temperature  conditions  than  the  winter  moisture. 
(Fletcher  and  Trlica  1980).  Regrowth  has  been  observed  after  late  season 
rainfall  (Cook  1977).  Growth  usually  begins  in  the  latter  part  of  March 
or  in  early  April  provided  soil,  moisture  and  weather  conditions  are 
favorable.  In  dry  summers,  plant  growth  stops,  the  herbage  becomes 
dry  and  brittle;  and  frequently,  the  plants  fail  to  produce  flowers  and 
seed.  (Hutchings  and  Stewart  1953). 
Phenology  (Pearson  1979) 

Observations  of  Indian  Ricegrass  indicate  that  the  plant  goes 
through  four  distinct  phenological  stages: 

1.  Sprouting  in  spring  as  dormancy  is  broken 

2.  Vegetative  growth  stage 

3.  Anthesis  and  maturation 

4.  .Dormancy  stage 

Growth  begins  when  the  soil  temperatures  at  the  15  cm  depth  warm  up 
to  4°C  and  remain  for  at  least  3  or  4  days.  During  the  sprouting  stage, 
foliage  production  appears  to  be  essentially  linear.  Maximum  plant 
height  was  also  correlated  with  soil  temperature. 

Irrigation  increased  the  amount  of  foliage  produced  and  delayed  the 
time  of  reaching  the  maximum  leaf  length. 

Stipa  Comata— Needle  and  Thread  Grass 

Stipa  comata  was  to  have  been  the  fourth  species  studied  but 
adequate  information  was  not  available  either  from  the  literature  or 
from  the  BLM  plant  surveys.  Crested  wheat  grass  was,  therefore, 
substituted.  By  the  time  the  decision  to  change  species  was  made, 
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there  was  not  time  to  make  as  complete  a  literature  survey  as  we  would 
have  liked. 

Agropyron  cristatum  or  desetorium — Crested  Wheat  Grass 

Crested  wheat  grass  begins  growth  early  in  the  spring,  but  the 
wheat  grasses  are  not  as  palitable  to  sheep  as  Bromus.  It  is  one  of  the 
first  grasses  to  start  growth  in  the  spring  and  continues  to  furnish 
pasture  until  late  in  the  fall  when  moisture  conditions  are  favorable. 
Ordinarily  the  poorest  growth  is  made  during  the  hot,  dry  periods  in 
midsummer.  For  revegetation,  the  seed  should  be  sown  as  early  in  the 
spring  as  possible  (Wilson  1931) 
Germination 

Good  germination  occurred  with  temperatures  of  20°C  and  30°C,  but 
the  highest  rate  of  germination  was  20°C.  The  seeds  can  germinate  under 
fairly  high  moisture  stress  if  the  temperatures  are  optional.  McGinniess(1960) 
Phenology 

Dr.  Fred  Provenza  (1980)  of  the  USU  range  department  who  is  working 
with  crested  wheat  grass  west  of  Eureka,  Utah  reports  that  crested  wheat 
grass  developed  into  the  boot  stage  in  his  study  area  during  1980  by  the 
3rd  week  in  May,  and  seeds  were  developing  well  by  the  first  of  June. 
The  plants  grew  quickly  this  year  during  the  summer. 

Production 

The  best  combination  of  temperature  and  precipitation  for  predicting 
unfertilized  mature  yield  was  July-May  precipitation  plus  mean  March,  April 
and  May  temperatures  at  the  Oregon  site  of  this  study.  (Sneva  1977)  The 
February  temperature  combined  with  March  precipitation  accounted  for  83% 
of  the  variability  in  spring  yield  at  this  study  site.  By  the  time  the 
mature  stage  was  reached  in  July,  additional  moisture  did  little  to  increase 
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the  yield. 

PROCEDURES  AND  TECHNIQUES 

As  was  indicated  previously,  the  work  on  modeling  of  tomatoes  and  the 
competing  weed  species  indicated  strongly  that  the  growth  curve  defined  by 
a  combination  of  equations  1  and  2  increased  the  predictive  accuracy 
considerably  over  the  results  achieved  with  equation  1  only.  Use  of  the  two 
equations  requires  the  determination  or  estimation  of  three  cardinal 
temperatures  for  each  species  rather  than  only  two. 
Cardinal  Temperature  Selection 

The  review  of  the  literature  previously  reported  on  failed  to  supply 
specific  cardinal  temperature  information  for  all  of  the  species  in 
particular  values  of  TC.  There  were  even  contradictions  on  the  values  for 
the  other  two  temperatures  but  it  appears  that  a  good  starting  point  for 
the  four  species  under  consideration  will  be  the  following:  TB  =  4°C  (39°F). 
TU  =  25°C  (77°F)  and  TC  =  36°C  (96.8°F).  The  36°C  value  was  selected  since 
in  all  of  the  modeling  of  this  type  to  date  the  rate  of  decrease  of  plant 
development  as  temperatures  increase  above  the  optimum  temperature  has  always 
been  much  more  rapid  than  the  rate  of  increase  in  development  as 
temperatures  increase  below  the  optimum  temperature. 

After  our  literature  review  and  finding  very  little  information  on 

the  cardinal  temperatures  for  the  different  plant  species,  we  Dropose  to  use 

the  following  general  categories  until  more  specific  values  can  be  measured: 

TABLE  2 
Cardinal  Temperatures  for  Three  Plant  Categories 
Plant  Type  TB  TU  TC 

Cool  Temperature  Plants  0°C  (32°F)  20°C  (70°F)  32°C  (90°F) 
Intermediate  Temp.  Plants  4°C  (39°F)  25°C  (77°F)  36°C  (97°F) 
Warm  Temperature  Plants     10°C  (50°F)   30°C  (86°F)   42°C  (108°F) 
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In  the  reviewed  literature,  several  of  the  plant  species  had  reports 
of  different  cardinal  temperature  classifications,  but  the.  strongest 
evidence  points  to  the  four  species  under  study  falling  in  the  intermediate 
category. 
Soil  Moisture  Content 

The  values  of  the  moisture  content  of  the  soil  at  the  wilting  point 
(WP)  and  field  capacity  (FC)  were  extracted  from  the  SCS  soil  survey  reports. 
The  values  for  the  soil  at  Price  were  fairly  site  specific,  but  data  of 
equivalent  accuracy  could  not  be  obtained  for  the  Vernal  site  since  the 
phenology  data  were  taken  over  a  wide  range  of  elevations  and  locations  in 
the  Vernal  area.  On  the  basis  of  the  general  sorl  classifications  in  the 
area,  it  appeared  that  the  same  values  could  be  used  for  both  sites  in  the 
preliminary  modeling  efforts.  Values  of  WP  =  55  mm  (2.17  inches)  and 
FC  =  127  mm  (5.0  inches)  were  selected.  Since  the  winter  of  1977-78  and 
1978-79  at  both  Price  and  Vernal  were  relatively  dry  it  was  assumed  that  the 
moisture  content  of  the  soil  at  the  beginning  of  the  growing  season  was 
only  50%  of  field  capacity  at  both  sites. 

As  constants  for  calculating  the  influence  of  moisture  on  the 
accumulation  of  GDH,  it  was  assumed  that  the  date  the  snow  cover  disappeared 
was  March  1,  (day  number  182  on  the  modified  Julian  calendar  which  begins 
September  1).  It  was  further  assumed  that  the  length  of  the  growing  season 
at  Price  was  140  days  and  at  Vernal  120  days  when  calculating  the  moisture 
usage  by  the  plants. 

There  was  not  sufficient  information  available  to  test  our  current 
method  of  calculating  evapotranspiration  and  the  influence  of  soil  moisture 
content  on  the  specific  species.  Surface  texture  characteristics  have  a 
significant  affect  on  available  soil  moisture  in  low  precipitation  climate 


16. 

regimes  (less  than  12  inches  annually).  The  General  soil  classification 
does  not  reflect  these  surface  texture  influences.  Some  actual  values  of 
the  amount  of  moisture  in  the  soil  are  required  to  refine  the  soil 
moisture  calculations.  These  improvements  may  modify  the  numerical 
value  of  the  GDH  constants  for  the  various  phenological  stages,  the  height 
and  the  biomass  calculations  developed  in  the  current  research;  but  once  the 
more  accurate  values  are  determined,  the  results  will  have  a  much  broader 
spatial  range  of  application. 
Estimating  the  Beginning  of  Growth 

No  information  on  the  time  growth  began  was  available  at  any  of  the 
observational  sites  for  any  plant  species.  It  was  therefore  necessary 
to  estimate  the  beginning  dates.  Winter  temperatures  in  most  areas  of 
Utah  are  too  cold  for  any  form  of  plant  growth.  The  assumption  was 
therefore  made  that  growth  began  on  the  date  in  the  spring  when  the  first 
accumulation  of  growing  degree  hours  for  the  season  was  accumulated  and 
followed -by  at  least  four  additional  days  with  temperatures  warm  enough  to 
continue  the  accumulation.  The  dates  of  beginning  of  growth  thus 
determined  are  shown  in  the  phenology  tables  in  the  appendix. 

Since  the  literature  review  indicated  that,  in  general,  the  species 
selected  for  study  all  would,  fall  in  the  intermediate  temperature  category 
of  TB  =  4°C,  TU  =  25°C  and  TC  =  36°C,  this  set  of  cardinal  temperatures 
was  used  in  the  modeling  work  being  reported. 
Calculating  GDH  data  and  Relating  to  Height,  Phenology  and  Production 

The  GDH  calculations  were  made  using  equations  (1)  and  (2)  and  the 
constants  discussed  previously.  A  sample  of  the  output  from  the  computer 
has  been  included  in  the  appendix  in  Table  23. 

To  develop  the  GDH  requirements  for  the  desired  phenological  dates, 
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the  accumulation  from  the  calculated  beginning  date  of  growth  to  each 
desired  stage  of  development  was  determined  for  each  available  year  and 
site.  The  accumulations  for  the  various  species  have  been  included  in 
the  appendix  tables  as  indicated.  The  requirements  for  the  various 
heights  and  various  determinations  of  biomass  were  determined  in  a 
similar  manner  and  have  also  been  included  in  the  tables. 

CONFIDENCE  LEVELS  ACHIEVED  FOR  OUTPUTS  BY  SPECIES 

The  data  that  was  available  from  the  BLM  areas  for  the  selected 
species  was  not  sufficient  to  run  any  valid  statistical  tests  to  -evaluate 
either  the  correlation  or  levels  of  confidence.  Graphs  which  show  the 
scatter  of  the  data  points  about  an  estimated  average  curve  have  been 
included  for  each  of  the  key  species  where  sufficient  information  was 
available. 

Discussion  of  Results  by  Species 
Stipa  comata— needle  and  thread  grass 

As  indicated  in  table  3A  in  the  appendix,  information  for  only  two 
years  was  available  for  this  species:  Vernal  1979  and  Vernal  1980. 
To  further  complicate  the  situation',  the  data  for  1979  was.  measured  at 
an  elevation  of  8600  feet  above  sea  level  while  the  1980  information  was 
observed  at  only  6000  feet  above  sea  level.  The  weather  station  at 
Vernal  is  located  at  an  elevation  of  5280  feet  above  sea  level. 

Orizopsis  hymenoides — Indian  Rice  Grass 

There  was  little  more  information  available  on  Orizopsis  hymenoides 
than  for  Stipa  comata.  Although  the  information  for  the  two  years  was 
collected  on  opposite  sides  of  the  basin,  they  were  observed  at  about  the 
same  elevation.  Only  three  height  measurements  were  taken  and  the  observer 
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indicates  that  she  is  quite  doubtful  of  the  validity  of  them. 

The  estimated  GDH  accumulations  for  the  various  phenol ogical  stages 
and  the  limited  height  and  biomass  information  have  been  tabulated  in 
Table  5,  6  and  7.  Again,  these  values  are  very  preliminary  due  to  the 
quality  and  limited  amount  of  basic  data.  A  very  preliminary  dryweight 
Biomass  versus  accumulated  growing  degree  hours  curve  has  been  plotted 
in  Figure  3.  There  is  considerable  scatter  to  the  few  points  but  the 
trend  is  readily  demonstrated. 

Agropyron  cristatum — crested  wheat  grass 

More  data  was  available  for  crested  wheat  grass  than  for  any  of  the 
other  species,  but.  again  there  is  a  problem  with  the  1980  information  at 
Vernal  being  taken  at  an  elevation  of  6600  feet  above  sea-  level  while  the 
observational  station  was  located  at  only  5280  feet.  The  basic  information 
has  been  tabulated  in  Table  3CX  and  3C2  in  the  appendix.  A  comparison  of 
the  accumulations  of  growing  degree  hours  for  the  various  phenological 
stages  indicates  that  the  coefficient  of  variation  is  reasonably  small, 
although  with  the  small  quantity  of  data  these  should  be  interpreted 
with  caution. 

The  coefficient  of  variation  (CV)  for  the  boot  stage  is  7.4%  but  with 
only  three  observations.  The  peak  flower  CV  =  2.0%  with  four  observations, 
the  seed  ripe  stage  CV  =  18.0%  with  4  observations,  and  for  the  seed 
disseminate  stage  CV  =  7.2%  again  with  4  observations. 

The  beginning  of  vegetative  growth  shows  a  large  difference  between 
the  accumulations  at  Price  and  the  accumulations  at  Vernal.  This  is  due 
to  a  difference  in  the  interpretation  of  the  definition  of  that  stage  of 
development  by  the  different  observers. 

The  fit  of  both  the  height  and  biomass  data  (Figures  4  and  5)  is  much 
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better  than  any  of  the  other  species,  but  there  is  still  considerable 
scatter.  However,  when  one  considers  the  problems  of  obtaining  reliable 
data  of  this  type,  the  fit  is  surprisingly  good. 

The  analysis  of  Agropyron  cristatum  indicates  the  potential  for 
application  of  this  form  of  modeling  to  the  determination  of  range 
development  and  production  quite  well.  If  we  eliminate  from  the  data 
those  points  which  show  a  decrease  in  production  after  a  larger  GDK 
accumulation,  the  correlation  coefficient  between  biomass  and  accumulated 
GDH  is  88%  with  a  standard  error  of  2209  GDH.  If  all  of  the  available 
data  points  are  considered,  the  correlation  coefficient  is  63%  with  a 
standard  error  of  8203  GDH.  Once  again  it  should  be  recognized  that 
such  statistics  are  not  statistically  conclusive  when  one  considers  the 
limited  amount  upon  which  they  are  based. 

Hilaria  jamesii— -gall eta  grass 

Data  on  galleta  grass  was  available  for  only  three  years:  1978  at 
Price,  1979  at  Vernal  and  1*980  at  Vernal.  Again,  we  were  unable  to 
utilize  1980  data  since  it  was  observed  at  an  elevation  of  only  4800  feet 
above  sea  level.  Due  to  the  nature  of  the  basin,  it  was  again  assumed 
that  the  Vernal  1979  data  would  be  fairly  compatible  even  though  the 
observations  were  taken  at  6,000  feet  above  sea  level.  In  determining 
the  GDH  requirements  for  the  various  phenological  stages,  we  again 
encountered  some  problems  in  the  vegetative  growth  stage  (see  Table  3Ul 
in  the  appendix).  There  was  also  a  difference  in  the  accumulated  GDH 
for  the  boot  stage;  but  the  peak  flowering,  seed  ripe  and  seed  disseminate 
accumulations  are  quite  close  for  the  Price  1978  and  Vernal  1979 
observations.  The  averages  required  for  the  various  stages  as  used  in  the 
predictive  equations  have  been  tabulated  in  Tables  5,  6  and  7. 
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The  graphs  of  Dry  weight  versus  GDH  and  Height  versus  GDH  are  shown 
in  Figures  6  and  7.  Again,  there  was  not  sufficient  data  for  valid 
statistical  tests;  but  the  degree  of  scatter  indicates  a  fairly  reliable 
relationship.  The  Biomass  data  for  1980  were  also  plotted,  but  it  appears 
that  the  information  as  tabulated  may  have  been  off  by  a  factor  of  10. 
The  biomass  data  at  Price  was  the  sum  of  the  total  weight  sampled  and  not 
that  of  individual  plants;  but  when  we  assumed  a  sample  of  10  plants,  the 
values  of  the  production  trends  were  similar.  Analyzing  the  five  samples 
during  1979  showed  a  correlation  coefficient  of  97.7%  and  a  standard 
error  of  566  for  this  year  which  was  the  only  one  with  what  we  considered 
as  reliable  data.  These  are  poor  statistics  but  does  give  a  rough  idea 
of  the  relationship  and  the  potential  for  predicting  biomass  of  this 
species. 

Atriplex  canescens--four  wing  salt  bush 

Only  one  year  of  information  was  available  for  this  species:  Vernal 
T979.  There  was  no  indication  of  the  elevation  at  which  these  observations 
were  taken,  but  the  information  was  evaluated  assuming  that  the  data  may 
have  been  representative  of  the  weather  station  elevation.  The  basic  data 
tabulation  can  be  found  in  the  appendix  Table  3E.  The  limited  Biomass 
and  height  information  has  been  plotted  in  Figures  8  and  9  .  The 
information  does  show  a  fair  fit  to  the  general  production  curve  found  in 
the  other  species.  This  is  of  course  only  a  'jery   rough  approximation  of 
the  actual  modeling  curves. 

Artemesia  Tridentata— Sage  brush 

Since  only  limited  amounts  of  data  were  available  for  the  selected 
species,  a  plot  was  also  made  of  sage  brush  since  the  elevation  of  the 
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observational  sites  for  this  species  were  fairly  representative  of  the 
weather  station  at  Vernal.  Here  again,  we  found  that  the  information  was 
not  very  satisfactory.  I  suspect  that  the  leader  length  data  at  Vernal 
during  1979  may  have  been  expressed  in  millimeters  instead  of  centimeters, 
but  I  can't  be  certain.  The  biomass  information  was  not  compatible  with 
the  various  sites "either;  hence,  the  plots  were  not  included  in  the 
summary.  The  phenology  data  appears  to  be  a  little  better  and  the  average 
GDH  accumulations  can  be  found  in  Table  9. 

COMPARISON  OF  OBSERVED  AND  CALCULATED  PRODUCTION  WITH  NORMAL  VALUES 

One  question  that  is  a  very  important  consideration  in  range 
management  is  "How  does  the  current  year's  production  compare  with 
normal  production  in  any  particular  area.  To  evaluate  the  use  of  the  model 
in  answering  this  type  of  question,  normal  temperature  and  precipitation 
values  for  both  Vernal  and  Price  were  entered  into  the  model  to  predict 
the  normal  production.  Two  runs  were  made.  The  first  assumed  that  the 
soil  moisture  at  the  beginning  of  the  growing  season  was  50%  of  field 
capacity,  the  same  values  assumed  in  developing  the  models.  The  second 
assumed  that  the  soil  moisture  was  100%  of  field  capacity.  The  50% 
run  resulted  in  predicted  production  not  too  different  from  the  values 
obtained  by  inputing  the  temperature  and  precipitation  for  the  specific 
years  for  which  data  was  available.  (Figure  10) 

When  the  observed  production  values  were  averaged  over  the  growing 
season  and  compared  with  predicted  normal  values  based  on  100%  of  field 
capacity,  the  1979  average  production  at  Price  was  83%  of  normal.  The 
1980  production  at  Vernal  was  88%  of  the  calculated  normal  production  at 
that  site.  The  1979  production  at  Vernal  was  107%  of  the  calculated 
normal . 
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In  determining  the  actual  production  figures,  those  observed 
values  which  were  followed  by  lower  production  measurements  were 
assumed  in  error  and  discarded  in  making  the  average  calculations  for 
the  season.  The  production  figures  at  Vernal  in  1979,  for  example, 
showed  much  larger  variations  than  either  of  the  other  season's  measurements. 
The  measured  values  have  been  tabulated  in  Table  17. 

A  second  test  compared  the  production  values  as  predicted  by  the 
models  for  each  year  with  the  predicted  normal  values  in  each  study  site. 
The  values  for  Agropyron  cristatum  and  Hilaria  jamesii  have  been  .plotted 
in  Figures  11  and  12.  It  will  be  noted  that  the  average  production 
as  calculated  for  each  year  is  less  than  the  calculated  normal  production. 
These  tests  indicate  that  with  more  accurate  dry  weight  biomass  information 
such  predicted  comparisons  will  be  very  practical . 

SUGGESTED  PROCEDURES  FOR  EXPANDING  PHENOCLIMATOGRAPHY  MODELS  TO  ALL 
MAJOR  OR  IMPORTANT  RANGE  LAND  SPECIES 

From  an  analysis  of  the  species  studied  in  the  current  project 
combined  with  other  modeling  efforts  during  the  past  few  years,  it  is 
now  possible  to  pinpoint  the  more  important  soil  and  atmospheric  factors 
which  relate  to  phenological ,  growth  and  production  of  specific  range 
species.  To  expand  the  modeling  program  and  increase  its  accuracy,  we 
need  to  collect  the  following  types  of  information  for  the  desired 
sites  and  species:  &**«*u^o  *#J  l^^^^A-  ^ /..*■'■■<  «<* 

A.  Soils  Data  ^^    ^  ...  \J„M  &    :J »  ~*    **' 

l»^7'liC"l«  Wilting  point  (determined  from  a  soil  ana-lysis)  ©  /■SaX^r'* 
,  ,/  ■'"   .-  L2.  Field  Capacity  (also  determined  from  a  soils  analysis)  1  **"*>r 
>V>  3.  Soil  moisture  content  for  at  least  three  depth  intervals 
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J  ? 


/,,  t»M.e.  0  to  30,  30  to  60,  60  to  90  cms.  (0  to  1 ,  1  to  2,  and 
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2  to  3  feet).  This  information  needs  to  be  collected  at 
the  beginning  of  the  growing  season  and  at  about  3  week 
intervals  during  the  season. 

4.  Daily  maximum  and  minimum  soil  temperatures  at  the  10  cm. 
(4  inch)  depth  throughout  the  year. 

B.  Atmospheric  data 

1.  Daily  maximum  and  minimum  air  temperatures  as  measured  in 
a  standard  instrument  shelter  during  the  entire  year. 

2.  Daily  precipitation  measurements  during  the  entire  year. 

3.  Date  the  snow  disappears  if  possible  each  year. 

C.  Plant  information 

1.  Phenological  stages  for  desired  species  including: 

a.  First  observed  signs  of  growth. 

b.  Green  Up. 

c.  Observable  stages  of  vegetative  growth  such  as  the 
number  of  leaves  and  date  they  occur.  These  data 
will  vary  with  the  species  under  study. 

d.  Fruiting  stages  such  as:  the  beginning  of  development 
of  the  infloresence,  boot,  first  bloom,  peak  bloom, 
seed  ripe,  seed  disseminate,  etc. 

e.  Timing  and  amount  of  fall  regrowth  and  determination 
of  how  much  carries  over  into  the  following  spring. 

2.  Height  information  or  leader  length  at  regular  intervals 
and  record  of  date  on  which  each  observation  was  taken. 

3.  Dry  weight  biomass.  Average  per  culm  or  individual  leader, 
Should  be  collected  at  regular  intervals.  Weights  must 

be  carefully  checked  and  the  date  of  each  observation 
recorded. 

4.  Estimates  of  average  percent  of  plant  cover  on  each 
observational  date.  These  data  are  required  to  estimate 
changes  in  effective  evapotranspiration  on  the  range. 

D.  Non-field  determinations  from  the  literature  or  laboratory  and/ 
or  growth  chamber  studies  as  needed. 

1.  Cardinal  temperatures  for  each  species. 

2.  Unique  growth  and  development  characteristics. 
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3.  Day length  sensitivity 

4.  Influence  of  critical  stress  factors  such  as  atmospheric 
stress,  soil  moisture  stress,  disease  stress,  insect 
stress,  grazing  stress,  etc.  on  plant  growth  and 
development. 

With  the  above  information  as  input  data  to  the  preliminary  computer 
program  which  has  already  been  developed,  the  GDH-phenology,  GDH-Height 
and  GDH-dryweight  biomass  models  can  be  developed  for  each  desired 
species.  Once  the  model  constants  have  been  d-etermined  and  tested  over  a 
diverse  range  of  climate  and  soil  conditions,  the  models  can  be  used  to 
predict  the  desired  variables  wherever  the  necessary  soil  and  climate 
information  is  available. 

Once  the  models  have  been  perfected,  it  will  be  possible  to  utilize 
climate  information  from  long  term  weather  stations  to  study  trends  in 
development  and  relate  such  trends  to  cultural  practices  on  the  range  and 
evaluate  the  influence  of  such  stress  factors  as  drought,  cold  temperatures, 
excessive  moisture,  overgrazing,  etc. 

FEASIBILITY  AND  USE  OF  PHENOCLIMATOGRAPHY  MODELS  IN  RANGE  MANAGEMENT 

The  manager  of  western  rangelands  is  continually  faced  with  a 
multitude  of  conflicting  requirements.  (Pressure  groups,  many  of  whom 
are  lacking  in  a  knowledge  of  the  basic  facts,  are  continually  pressing 
for  utilization  of  the  ranges  in  a  manner  which  they  feel  will  best  suit 
their  own  personal  desires.  Many  such  pressures  are  quite  self-centered, 
selfish,  and  personally  motivated.) 

Most  of  the  current  decisions  which  must  be  made  involve  elements 
that  create  an  impact  on  the  air,  soil  and  water  aspects  of  the  environment 
and  may  have  yery   important  influences  on  the  welfare  of  plants  and  animals 
on  the  range.  (Such  factors  are  generally  interactive  and  present 
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knowlege  on  how  to  predict  such  interactions  is  far  from  adequate.) 

Perhaps  the  most  important  use  of  the  models  is  the  prediction  of 
phenological  development  and  production  for  use  in  grazing  management. 
Grazing  treatments  are  usually  keyed  on  plant  development  and/or  an 
allowable  utilization  of  current  year  production  of  key  species.  In 
the  past,  timing  of  grazing  treatments  has  been  designed  around  dates 
representing  plant  development  and  production  estimates  for  a  "normal" 
year.  To  achieve  maximum  plant  response  or  benefit  to  the  user,  grazing 
management  should  be  flexible  enough  to  allow  treatments  to  be  initiated 
as  closely  as  possible  to  the  actual  occurrence  of  the  phenological  stage 
or  utilization  level  desired. 

(Modeling  cannot,  of  course,  replace  actual  field  observations  by 
the  trained  range  manager;  but  it  can  provide  a  valuable  tool  for  increasing 
the  scope  and  efficiency  of  limited  field  observations  for  making  broader 
management  decisions.) 

Another  important  use  of  the  models  is  in  prediction  of  range 
readiness  and  production.  When  sufficient  climate/weather  information 
is  available  to  establish  long  term  averages  or  normal  temperature  and 
moisture  data,  the  current  year's  information  is  used  as  far  as  available 
and  the  assumption  is  made  that  beyond  this  period,  normal  temperatures 
and  moisture  conditions  will  occur.  Using  this  assumption,  it  is  possible 
to  predict  each  year  at  any  desired  time  the  stage  of  development  of  the 
desired  species  and  the  estimated  production.   (It  is  also  possible  to 
vary  the  estimates  in  terms  of  departures  from  normal,  either  above  or 
below,  during  the  period  when  normals  are  used  in  the  predictive  model. 
As  new  current  weather  information  becomes  available,  the  predictions 
can  be  updated  as  frequently  as  desired  to  obtain  more  recent  predictive 
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results  for  guiding  the  range  manager  in  his  decision  making.) 

A  third  application  of  the  models  is  to  use  them  as  a  form  of 
environmental  control  of  weather  factors  in  field  and  range  experiments. 
In  the  past  the  value  of  many  field  experiments  was  limited  because  the 
variable  nature  of  weather  and  its  influence  at  the  time  the  experiments 
were  performed  was  not  usually  evaluated.  Even  in  their  current  juvenile 
state  of  development,  these  models  will  enable  researchers  to  eliminate 
to  a  certain  extent  the  influence  of  weather  conditions  on  field  experiments 
by  enabling  researchers  to  compare,  conditions  at  the  time  of  their  own 
research  with  the  development  which  would  occur  under  normal  conditions 
or  in  any  other  particular  year  or  season.  Researchers  can  thus  define 
more  closely  the  impact  of  other  influences  which  might  be  modifying  the 
growth  and  development  of  the  plants  since  the  influence  of  variable 
weather  conditions  can  be  partially  eliminated. 

The  models  also  permit  the  range  manager  to  evaluate  the  influence 
on  a  longer  time  scale  of  various  cultural  and  management  decisions. 
For  example,  the  production  on  the  range  as  evaluated  by  the  model  is 
plotted  in  comparison  with  the  observed  information  over  a  period  of 
years.  It  is  then  possible  to  detect  trends  in  departures  from  the 
potential  production  or  development  which  can  in  many  cases  be  related 
to  certain  management  decisions. 

It  is  also  possible  by  use  of  the  models  to  calculate  production 
over  the  full  period  of  weather  records,  to  compare  these  with  available 
observed  production  data  and  to  analyze  trends  in  production  over  a 
longer  period  of  time  than  actual  production  data  may  be  able  to  supply. 
Such  a  study  will  enable  the  range  manager  to  better  understand  the 
potential  influence  on  production  of  future  periods  of  drought,  excessive 
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moisture,  abnormally  cold  or  warm  temperatures  or  current  grazing 
practices  on  total  range  production.  The  range  manager  can  further  be 
better  prepared  to  develop  contingency  plans  for  reducing  the  economic 
impact  of  such  extreme  weather  occurrences  on  the  users  of  the  range. 
The  immediate  effects  of  the  drought  or  other  abnormal  weather  occurrences 
can  also  be  quantified  and  analyzed  as  they  relate  to  adjustments  in 
livestock  numbers,  periods  of  use  or  supplemental  feed  requirements. 

Since  the  models  can  predict  the  influence  of  varying  temperatures 
and  moisture  regimes  on  range  production,  even  if  not  perfectly,  it  is 
also  possible  to  utilize  them  in  predicting  the  economic  values  of  such 
practices  as  addition  of  fertilizer,  reduction  of  competing  of  nondesirable 
plant  species  and  the  impact, of  increased  insect  populations  on  production. 
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TABLE  3 


DATA  AVAILABLE  FOR  MODEL 
DEVELOPMENT  AND  TESTING  FROM  PRICE,  UTAH 


Species 

Years 
1978 

Phenology 

Hei 
Leader 

ght  or 
Lengths 

Biomass 
Total 

*AGDE2 

X 

X 

'.  X 

*AGDE2 

1979 

X 

X 

• 

*HIJA 

1978 

X 

X 

X 

ARTR2 

1978 

X 

X 

X 

ARTR2 

1979 

X 

X 

OPPO 

1978 

X 

X 

OPPO 

1979 

X 

X 

SAKA 

1978 

X 

X 

SAKA 

1979 

X 

X 

*These  were  the  only  data  for  the  desired  key  species  to  be  modeled. 


TABLE  4 


DATA  AVAILABLE  FOR  MODEL 
DEVELOPMENT  AND  TESTING  FROM  VERNAL,  UTAH 
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Species 


Site 
Elevations*   Years 


Phenology 


HIJA 

5800 

1979 

X 

HIJA 

4800 

1980 

X 

ORHY 

5800 

1979 

X 

ORHY 

6000 

1980 

X 

SAKAT 

5100 

1979 

X 

AGCR 

7000 

1979 

X 

AGCR 

6600 

1980 

X 

STC04 

3600 

1979 

X 

STC04 

6000 

1980 

X 

Height  or 
Leader  Lengths 


Biomass 
Per  Culm 
Or  Leader 


X 
X 

X 
X 


X 
X 

X 
X 


ATCA 


5800? 


1979 


*The  weather  station  at  Vernal  is  located  at" an  elevation  of  5280  feet  above 
sea  level,  but  reading  from  the  topographic  maps  the  locations  given, 
the  data  was  collected  at  elevations  ranging  from  4800  feet  above  sea  level 
to  8600  feet;  and  in  some  cases,  on  different  sides  of  the  Uinta  Basin. 
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TABLE  5 


Grasses:  GDH  accumulation  for  related  phenological  stages 


Phenology  stages 

1. 

Begin  Growth 

2. 

Two  Leaves 

3. 

Three  Leaves 

4. 

Four  Leaves 

5. 

Five  Leaves 

6. 

Boot 

1. 

First  Flower 

8. 

Peak  Flower 

9. 

Peak  Seed  Ripe 

0. 

Seed  Dissemination 

ORHY 

AGCR 

HIJA 

0 

0 

0 

1100 

2200 

1900 

3700 

3100 

5200 

5000 

6900 

6425 

6309 

7006 

9188 

9571 

9927 

11622 

12802  . 

11716 

12675 

16933 

14402 

TABLE  6 


Grasses:  GDH  accumulation  for  related  heights 


Heights  (cm) 


10 


15 


20 


25 


30 


35 


40 


HIJA 
AGCR 
ORHY 


2600   5900   9300 

1100   2500   3700   5000   6700   9200 

1500   3600   7000 
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TABLE  7 


Grasses:  GDH  accumulation  for  selected  plants  for  dry  weight  biomass 


Biomass  (gms)    5  10  15  20 

AGCR  1300  2700  4300  6200 

HIJA*  4700  8300  11000  13000 

ORHY  2000  3500  4900  6000 


25    30    35 
8500  11400  14500 


40 


45 


50 


6900   7600   8200   8700   9300   9800. 


*  Biomass  X  10"  grams 


The  growth  of  the  seed  stalk  or  infloresence  apparently  does  not  begin 
until  the  grass  plant  reaches  a  certain  stage  of  development.  In  terms  of 
growing  degree  hours,  the  required -accumulation  before  the  plant  is  sensitive 
to  daylength  as  a  trigger  is  on  the  order  of  2000  to  2500  GDH.  This 
accumulation  seems  to  correspond  to  a  stage  of  development  ranging  between 
3  leaf  and  5  leaf  stages  depending  upon  the  grass  species. 

Utilizing  the  above  assumption,  the  following  table  shows  the  relationship 
between  th'e  height  of  the  infloresence  and  the  accumulation  of  GDH.  Development 
of  the  infloresence  is  assumed  as  0  em's  at  2500  GDH. 


TABLE  8 


Grasses:  GDH  accumulation  compared  to  height  of  infloresence 


Heights    5  10  15  20  25  30  35  40    45    50    55 

AGCR    2800  3100  3500  4100  5000  5900  6900  8100  9500  11200  13200 

HIJA   .3300  4300  5400  6400  7500  8700  10100  11700 

ORHY    3100  3800  4600  5800  7000  8300  10000  12200 
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TABLE  9 
Shrubs:    GDH  accumulation  for  selected  phenological  stages 
Phenological  stages 

1.  Begin  growth 

2.  Green  Up 

3.  Vegetative  Stage 

4.  Floral  Bud  Development 

5.  First  Flower 

6.  Peak  Flower 

7.  Seed  Ripe 

8.  Seed  Dissemination 


TABLE  10 


Shrubs:   GDH  accumulation  for  selected  plant  leader  length 


Lengths  (Cms) 

10 

20 

30 

40 

50 

60 

70 

ATCA 

800 

1600 

2500 

3400 

4700 

6100 

8100 

Based  on  very   limited  information.  The  lengths  reported  seem  to  be  out 
of  line  with  more  normal  plant  growth  for  the  species. 


TABLE  1 1 
Shrubs:   GDH  accumulation  for  selected  pi  ant- dry  weight  biomass 


Biomass  (gms) 

5 

10 

15 

20 

25 

30 

ATCA*  • 

2000 

3900 

5900 

8000 

10800 

14900 

These  weights  are  quite  questionable  but  do  follow- a  good  curve. 


TABLE  12 

ARTEMESIA  TRIDENTATA 
Phenology,  Height  and  Biomass 


Begin   Green   Vegetative   Floral  Buds   Start 
Growth   Up     Growth      Develop    Flower 


Peak     Peak       Seed 
Flower   Seed  Ripe   Disseminate   Dormancy   Elevation 


Vernal  '79  Dates 

3/13 

5/22 

ACC  GDH 

0 

5,295 

Dry  Weight 
Biomass/leader 

14  gms 

Avg  leader 
length 

82  cms 

Vernal  '80  Dates 

2/14 

4/29 

5/20 

ACC  GDH 

0 

3,228 

5,457 

Dry  Weight 
Biomass/Leader 

.29  gms 

.09  gms 

6/5 
7,500 
24  gms 

67  cms 


7/21 
14,350 
.30  gms 


9/25 
23,988 
18  gms 

61  cms 


10/18 
26,560 
29  gms 

79  cms 


5,800 


6,000 


Avg  leader  length 


Price  '78  Dates 

2/26 

.4/29 

5/22 

ACC  GDM 

0 

3,000 

5,352 

Price  '79  Dates 

2/28 

4/20 

ACC  GDH 

2,090 

7/11 
13,406 

10/9 
26,435 

10/25 
28,839 

11/7 
29,610 

7/1 
11,605 

9/25 
25,706 

Price  '78  Dates 

2/26 

5/22 

ACC  GDH 

0 

5,352 

Biomass 
Dry  Weight 

0.9  gms 

Dry  Weight  Biomass/Leader 
6/19       7/5     7/18 

9,746  12,405   14,490 

1 .4  gms  1.3  gms  0.9  gms 


7/25 

8/8 

8/22 

15,616 

17,662 

19,928 

0.9  gms 

1 . 5  gms 

1 . 4  gms 

CO 
CTv 


TABLE  13 


ORYZOPSIS  HYMENOIDES  (Indian  Rice  Grass) 
Phenology  and  Dry  Weight  Biomass 


Begin 
Growth 


Green 
Up 


Vegetative 
Growth 


Vernal  '79  Dates    3/13 
ACC  GDH 
Wt  of  Culm 
Height 

Vernal  '80  Dates   2/14 
ACC  GDH 
Wt  of  Culm 
Height 


Boot 


First 
Flower 


5/22 

5/30 

5,295 

6,612 

18  gms 

35  gms 

12  cm 

5/13- 

5/27 

4,765 

6,237 

7  gms 

21   gms 

Peak 
Flower 


Peak 
Seed  Ripe 


Seed 
Disseminate 


Dormancy      Elevation 


6/25  7/9 

10,108  12,211 

59  gms  80  gms 

18  cm  40  cm 

6/11  6/30 

8,267  11,034 

36  gms  30  gms 


5,800 


7/10 
12,675 


6,000 


CO 


TABLE  14 


STIPA  COMATA  (Needle  and  Thread  Grass) 
Phenology  and  Dry  Weight  Biomass 


Begin 
Growth 


Green 
Up 


Vegetative 
Growth 


Boot 


First 
Flower 


Peak 
Flower 


Peak 
Ripe 


Seed 
Disseminate 


Dormancy       Elevation 


Vernal  '79  Dates 

4/13 

6/6 

6/14 

ACC  GDH 

0 

7,500 

8,491 

Wt  of  Culm 

7  gms 

23 -gms 

Height 

24  cm 

Vernal  '80  Dates 

2/14 

4/29  • 

5/27 

ACC  GDH 

0 

3,228 

6,237 

Dry  Wt  of  Culm 

14  gms 

18  gms 

Height 

7/6  7/12 

11,790  12,644 

36  gms  27  gms 

18  cm  34  cm 

6/24  6/30 

10,144  11,034 

40  gms  40  gms 


7/27 
14,956 

30  gms 


7/10 
12,675 

38  gms 


8,600 


6,000 


8 


TABLE  15 

HILARIA  JAMESII  (Gall eta  Grass) 
Phenology  and  Dry  Wt  Biomass 


Begin 

Green 

Vegetative 

Start 

Peak 

Peak 

Seed 

Growth 

Up 

Growth 

Boot 

Flower 

Flower 

Seed  Ripe 

Disseminate 

Dormancy 

Elevation 

Price  '78  Dates 

2/26 

4/2 

5/18 

6/2 

6/19 

6/28 

7/18 

9/25 

6,000 

ACC  GDH 

0 

1,238 

4,906 

6,971 

9,746 

11,220 

14,490 

24,936 

Dry  Wt  Biomass 

Total 

5/25 
13. 5  gms 

16.2  gms 

16.4  gms 

16.5  gms 

Vernal    '79  Dates 

3/13 

5/30. 

6/18 

6/25 

7/9 

7/23 

5,800 

ACC  GDH 

0 

6,612 

9,106 

10,108 

12,211 

14,314 

Dry  Wt 
Biomass/Culm 

0.7  gms 

1 . 2  gms 

1 . 2  gms 

2.0  gms 

2.2  gms 

Vernal    '80  Dates 

2/14 

5/7 

5/14 

5/22 

6/12 

6/25 

4,800 

ACC  GDH 

0 

4,298 

4,850 

6,346 

9,253 

10,277 

Dry  Wt 
Biomass/Culm 

■ 

o06  gms 

. 05  gms 

.20  gms 

.16  gms 

CO 

to 


TABLE  16 
HILARIA  JAMES  1 1 
Height  Data 

Price 

1978 

Dates 

4/4 

4/25 

5/8 

5/18 

6/2 

6/19 

6/28 

7/18 

7/25 

8/8 

8/22 

9/12 

9/25 

10/25 

Leaf 
Stage 

1-2 
Leaves 

2-3 
Leaves 

3  leaf 

3  leaf 

5  leaf. 

5  leaf 

Flowers 

Beginning 

to  dry 

Seeds 
Mostly 

Gone 
6  leaf 

Leaves 

Drying 

Out 

30-35% 
Dry 

No  Growth 
Still 
Green 

Hear  Wash 

Varying 
Amounts  of 
Regrowth 
Near  Wash 

No  change 

Average  Ht 

1-2  cms 

4  cms 

5-7  cms 

7-10  cms 

inflores- 
cence 
Ht  20  cms 

Leaf  Ht 
15-17  cms 
inflores- 
cence 
30  cms 

inflores- 
cence 

stalks 

drying 
out 

15-17 
cms 
Ht. 

• 

ACC  GDH 
5800 

1,334 

3,000 

3,739 

4,906 

6,971 

9,746 

11,220 

14,490 

15,619 

17,662 

19,928 

23,480 

24,936 

28,839 

Vernal 

1979 

Dates 

5/30 

6/25 

7/9 

7/23 

Leaf  Stage 

Probably 
Ht  of 

inflores- 
cence 

Probably 
Ht  of 

inflores- 
cence 

Probably 
Ht  of 

inflores- 
cence 

Probably 
Ht  of 

inflores- 
cence 

Average  Ht 

18.3  cms 

24.4  cms 

45.7  cms 

15.2  cms 

ACC  GDH 

6,612 

10,108 

12,211 

14,314 

' 

TABLE  17 

AGROPYRON  CRISTATUM  (Crested  Wheat  Grass) 
Phenology  and  Dry  Weight  Biomass 


Begin 
Growth 

Green 
Up 

Vegetative 
Growth 

Boot 
Stage 

Start 
Flower 

Peak 
Flower 

Peak 
Seed  Ripe 

Seed 
Disseminate 

Dormancy 

Elevation 

Price  '78  Dates 

2/26 

3/25 

5/25 

6/2 

6/19 

6/28 

8/8 

9/25 

6,000 

ACC  GDH 

0 

539 

5,781 

6,971 

9,746 

11,220 

17,662 

24,936 

Dry  Weight 
Biomass 

17.8  gms 

20.7  gms 

26.8  gms 

* 

Price  '79  Dates 

3/3 

4/1 

4/5' 

6/5 

6/19 

6/26 

7/23 

8/22 

6,000 

ACC  GDH 

0 

368 

496 

7,661 

9,724 

10,812 

15,117 

19,901 

Dry  Weight 
Biomass 

—  — 



"*  — 

-- 

_.- 

— 

Vernal    '79  Dates 

3/13 

4/24 

5/23 

5/31 

6/20 

7/17 

8/13 

5,800 

ACC  GDH 

0 

2,532 

5,454 

6,682 

9,391 

13,347 

17,434 

. 

ACC  Biomass 
Dry  Weight/Culm 

18  gms 

35  gms 

22  gms 

35  gms 

35  gms 

41   gms 

Vernal    '80  Dates 

2/15 

4/29 

5/14 

5/29 

6/19 

7/31 

8/12 

6,600 

.ACC  GDH 

0 

3,229 

4,850 

6,463 

9,424 

15,829 

17,520 

ACC  Biomass 
Dry  Weight/Culm 

18  gms 

17  gms 

17  gms 

22  gms 

34  gms 

20  gms 

m 


TABLE  18 
AGROPYRON  CRISTATUM 
Height  Data 


.       , 

Price  '78 

Dates 

Begin 
2/26 

3/21 

4/4 

4/25 

5/8 

5/18 

6/2 

6/19 

6/28 

7/18 

7/25 

8/8 

8/22 

9/12 

9/25 

Leaf 
Stages 

2  leaves 

3  leaf 

4  leaf 

4-5  leaf 

5  leaf 

in 
flower 

full 
flower 

seeds 
develop 

peak 
seed 
ripe 

peak 
seed 
ripe 

seeds 
dissem- 
inate 

seeds  dis- 
seminate 
leaves  dry 

seeds  dis- 
seminate 
culm  green 

plants 
dry 

Average 
Ht. 

5  cms 

5-6  cms 

12  cms 

15-20  cms 

25-30  cms 

inflores- 
cence 
40  cms 

leaf  30cms 
inflores- 
cence 
40-45  cms 

inflores- 
cence 
45-50  cms 

inflores- 
cence 
48-54  cms 

inflores- 
cence 
60  cms 

ACC  GDH 

0 

228 

1334 

2908 

3739 

4906 

6971 

9746 

11,220 

14,490 

15,619 

17,661 

19,928 

23,480 

24,936 

Price  '79 
Dates     '■ 

3/3 

4/1 

4/26 

5/21 

6/19 

7/3 

7/19 

7/31 

8/22 

9/10 

Leaf     ■ 

Stages 

3-4  leaf 

3-4  leaf 

full 
flower 

seeds 

develop 

seeds 
dissem- 
inate 

seeds 
dissem- 
inate 

Average 
Ht 

10-15-  cms 

15-25  cms 

5-61  cms 

ACC  GDH 

0 

369 

2909 

5498 

9725   • 

11,958 

14,436 

16,334 

19,901 

23,126 

Vernal  '79 

Dates 

4/24 

5/23 

5/31 

6/20 

7/17 

8/13 

Leaf 
Stages 

Boot 
Stage 

Peak 
Flower 

Peak 
Seed 
Ripe 

Seed 
Dissem- 
inate 

Average 
lit 

12  cms 

46  cms 

70  cms 

55  cms 

ACC  GDH 

2532 

5454 

6682 

9391 

13347 

17434 

TABLE  19 


ATRIPLEX  CAhESCENS  (Four  Wing  Salt  Bush) 
Phenology 


Begin 
Growth 


Green 
Up 


Vegetative 
Growth 


Floral   Buds 
Develop 


Peak 
Flower 


Seed 
Ripe 


Seed 
Disseminate 


Vernal   '79  Dates 

2/13 

6/5 

7/15 

9/5 

9/25 

ACC  GDH 

0 

7,343 

11,621 

21,053 

23,988 

Dry  Wt/leader 

17  gms 

"  28  gms 

28  gms 

33  gms 

Length  of  leader 

67  cms 

79  cms 

46  cms 

10/12 
26,080 
35  gms 


TABLE  20 


Species  Name: 

SEED 

Description: 


GERMINATION  REQUIREMENTS 


Scarification 


Vernalization  (chilling) 


Soil  Temperature 


Soil  Moisture 


Light  Intensity 


Light  Duration  (day  length) 


Light  Quality 


Other 


Annual 


,  Perenial,  Grass,  Forb,  Shrub,  Tree  (Underline) 


Lower  Lethal 


Base 


Optimum 


Critical 


Upper  Lethal 


Notes: 


TABLE  21 


Species  Name: — — 

VEGETATIVE  STAGES 

Name  and  description  of  recognizable  stages  of  development 


Climate  Factors 

Lower  Lethal 

Base 

Optimum 

Critical 

Upper  Lethal 

■ 

t  iaht  Intensity 

. 

1 iaht  Oual itv 

■ 

1 . 

.—..——- 

Nutrient  Requirements 

__M. 

Other 

Other 

' 

-^ 

Notes  on  back 


TABLE  22 


s~\ 


Species  Name: 


FRUITING  STAGES 

Names  and  descriptions  of  recognizable  stages  of  development; 


Climate  Requirements 

Lower  Lethal 

Base 

Optimum 

Critical 

Upper  Lethal 

and  triggers 
Temperature 

Ho  situ re 

Humidity 

Daylength 

Light  Intensity 

Light  Quality 

i 

Other 

Notes: 


en 


47 


TABLE  23 

SAMPLE  OF  COMPUTER  OUTPUT 

SPECIES:  STCO 

DAY  #  OF  PHENOL.  STAGE  DEVEL. 
Stage  #  12345678 

Predicted  Day       0      0    260   282      0    308    314    325 

DAY  #  OF  HEIGHT  LEVELS  (CM) 
Height  5    '10    15    20     25     30     35     40 

Predicted  Day      00000000 

> 

DAY  #  OF  BIOMASS  DEVEL.  (GMS) 
Dry  weight 
Biomass  5      10    15    20     25     30     35     40 

Predicted  Day     247    265    275   284    292     298    303    308 


TABLE  24 
Comparison  of  Predicted  Actual  and  Normal  Dates  For 
AGROPYRON  CRISTATUM 


Phenology 


Begin 
Growth 

Leaf 

Leaf 

Leaf 

Leaf 

Boot 

First 
Flower 

Peak 
Flower 

Peak 
Seed  Ripe 

Seed 
Disperse 

Vernal  1977-78 
Vernal  1978-79 
Vernal  Normals 

AW  -  50% 
AW  -  100% 

216 

234 

232 
221 

227 

244 

244 
231 

249 
265 

257 
244 

265 
271 

272 
262 

276 
282 

282 
273 

302 
304 

304 
295 

325 
326 

325 
316 

351 

353 

351 

343 

Price  1978-79 

Price  Normals 
AW  =  100% 

233 
215 

242 
226 

258 
240 

Height 

269 
257 

of  Culm 

279 
267 

301 
288 

321 
307  , 

347 
332 

Vernal  1977-78 
Vernal  1978-79 

;. 

5  cm 
216 
234 

10  cm 
237 
254 

15  cm 
255 
263 

20  cm 
265 
271 

25  cm 
279 
285 

30  cm 
300 
302 

Vernal  Normal 
AW  =  50% 
AW  =  100% 

232 
221 

251 
237 

262 
250 

272 
262 

285 
276 

301 
293 

Price  1977=78 

Price  Normal 
AW  =  100% 

233 
215 

252 
233 

261 

246 

269 
257 

282 
269 

298  • 
285 

Vernal  1977-78 
Vernal  1978-79 
Vernal  Normal 
AW  =  100% 

Price  1978-79 
Price  Normal 

AW  =  100% 


TABLE  25 
Comparison  of  Predicted  Actual  and  Normal  Dates  For 

HILARIA  JAMESII 

Phenology 


Begin 

First 

Peak 

Peak 

Seed 

Growth 

Leaf 

Leaf 

Leaf 

Leaf 

Boot 

Flower 

Flower 

Seed  Ripe 

Disperse 

Vernal  1977-78 

233 

255 

267 

280 

281 

304 

316 

335 

Vernal  1978-79 

194 

247 

263 

272 

286 

287 

307 

319 

336 

Vernal  Normals 

AW  ■  100% 

234 

250 

264 

277 

278 

298 

309 

327 

Price  1978-79 

245 

261 

270 

283 

284 

303 

314 

331 

Price  Normals 

AW  =  100% 

230 

246 

258 

271 

272 

290 

300 

316 

0.5  cm 

1.0  cm 

1 .5  cm 

239 

274 

300 

255 

278 

303 

239 

270 

294 

254 

276 

299 

234 

264 

286 

Height  of  culm 

2.0  cm        2.5  cm        3.0  cm 


KO 


Figure  3 
ORIZOPSIS  HYMENOIDES 


50 


65 

60 

55 

50 

45 

40 

— 

35 

30 

25 

- 

20 

- 

15 

- 

10 

- 

5 

- 

s 

.-VERNAL,  UT  1979     5,800  feet  elevation 
x-VERNAL.  UT  1980     6,000  feet  elevation 


/. 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


/ 


s 


/ 


s 


1    I    I I I L 


J L 


10    20    30    40    50    60    70    80    90   100   110   120    130 
ACCUMULATED  GROWING  OEGREE  HOURS  X  10 

Figure  4 
A6R0PYR0N  CRISTATUM  (Crested  Wheat  Grass)  AGCR 
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Figure  5 
A6R0PYR0N  CRISTATUM  (Crested  Wheat  Grass)  AGCR 
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Figure  6 
HILARIA  JAMES I I  (Galleta  Grass)  HIJA 
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Figure  7 
HILARIA  JAMESII    (Galleta  Grass)   HIJA 
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Figure  8  . 
ATRIPLEX  CANESCENS   (Four  Wing  Salt  Bush)  ATCA 
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Figure     9 
ATRIPLEXCANESCENS  (Four  Wing  Salt  Bush)  ATCA 
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Figure  10 
AGROPYRON  CRI-STATUM  (Crested  Wheat  Grass) 
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Figure  11 

AGROPYRON  CRISTATUM 

DRY  WEIGHT  BIOMASS  vs   PREDICTED  DATES 
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Figure  12 

HILARIE  JAMESII 

DRY  WEIGHT  BIOMASS  vs   PREDICTED  DATES 


i 


2.25 
2.0 

1.75 

1.50 

1.25 

1.0 

0.75 

0.5 

0.25 


.  Vernal  (1977-7U) 
®  Vernal  (H7S-79) 
x   Price  (1973-73) 

CD  Vernal  Normal 
+  Price  Nomial 


CO  f   -o 


1 


I 


220   230   240   250   260   270   230   290   300   310   320   330   340   350 
OAY  NUMBER  IN  FRUIT  TREE  OAYS 


XORKFILEl  RANGE  (12/10/80) 


ll«6  PH  THURSDAY,  OECEMBER  11,  1980 


100 
200 
300 
100 
S00 
600 
TOO 
750 
600 
")00 
1000 
UOO 
1200 
1300 
1100 
1500 
IfaOO 
1700 
1600 
19  00 
200O 
2100 
2200 
2J00 
2100 
2500 
2bO0. 
270O 
2800 
2900 
3000 
3100 
3200 
3300 
3100 
3500 
3600 
3700 
3600 
3^00 
1000 
1100 
1200 
<l300 
1100 
150  0 
5500 
56iio 
5  7  00 
5000 
5  90  0 
6000 
6100 
6200 
6300 
6100 
6500 
M.QQ 


C  — 

c— 
c— 
c— - 

c— 

FILE 
FILE 
FILE 
FILE 
FILE 


THIS  IS  A  PROGRAM  OF  A  RANGE  MODEL  USED  TO  PREDICT  PHENOLOGY. 
GRO-TH  AMD  PEAK  ABOVE  GrIOUNO  SEASONAL  PRODUCTION 


5(KINDoREH0TE) 

6(KInD«REM0TE) 

T(KIM£»cPRINTE") 

8{KI'iD:PACH,FlLETYPEs8iTITLE»'VRN078') 

l3(Mt,-D  =  PACK,FILETYPE«8,TnLEJ,,VrtNCAL'> 
DIMENSION    LAN(2),NC0N(2),NAN<2),LCQN(2),IXNTC«75i2)»IXNP{«7S) 

COMMON  /PLK1/  XNP(<I75),XNT(175,2),ETTC«7S),EAU75) 

COHHon  /BLK2/  NG(175),SPK(2/10),NDG0H(175) 

COMMON  /DLM/  AH(«7S),TSAT(1) 

20  WRITE<6,10) 

io  for-utp    ','  please  enter  the  station*'/) 

REAO(5,700)    1ST  , 

700  FORliAT(Ifc) 

IF(IST.EO.0>G0    TO    S-m 
DO    29    M=1,U75 
ETT(M)s0.O 
E A  Ci  1  =  0.0 
29      CONTINUE 

*H1TEU,850) 

850  FORMATC  V  ENTER  B£G.  YEAR,  IF  IT  IS  A  LEAP  YCAR{YES»1', 
*i  ,0=0),'/'  DATA  TYPE  (HETRICol,  NONHETRIC»0) </) 

READ(5,851)  1Y0,1LEAP,IDTYPE 

851  F0RHAT(I«»2Il) 

c— -- 

C— .  OPTION  FOR  DAILY  OUTPUT  OR  PHENO.  STAGE  PREDICTION 

WRITE (6,26) 

26  FORMAT! '  ','  ENTER  OPTION  FOR  DAILY  OUTPUT  OR  PHENO.1, 

»l  PREDICTION  TABLESi/i  (DAILY  OUTPUT  ONLY«l , PHENO  OUTPUT  OnLY«0,«, 
*'  OR  C0THa2)</) 

READ(5,7I0)  IDOP 

FORMAT  (ID 


710 

C 

C 

c 


INPUT  SPECIE  AND  CARDINAL  TEMPS, 


WR1TE(6,25)  ,   „ 

25   ForMatI'  ','  ENTER  SPECIE  NAHE  AND  NUMBER,  ALSO  CARDINAL1* 
»t  TEMPS.  IN  A  A6,1I2  FORMAT!/) 
READ(5,27)  SPP,ISP,TL,TO,TU 

27  F0RH»Y(Ab,l2,JF210) 

IF  (ISP.E:).0)  CO  TO  9999 
23   wHITr.{6,B51)  .  • 

851  F0a"ATC  V  ENTER  CONSTANT  FOR  CALC,  OF  STRESS  TACTOR'/) 
READ<5,855)  C^NST 

655  F0fiMATtF5.il 
MITE  (6,652) 

852  FOSMATC  ','  IS  SOIL  MOISTURE 
*/'  CALCULATED  AND  AP?LIED(ANS, 

REAi)(5,851)  N0S.1,N0ATH 
8SJ  F0RHATI2II) 
NDT:156 
IF  (lLEAP.EO.I)  NDT*157 


OR  ATMOSPHERIC  STRESS  TO  6E«, 
EACH  dlTH  YESal,  NO=0)«/J 


00000100 

U0000200 

00000300 

000C0100 

00000500 

OOOOObOO 

0000*700 

00000750 

OOOOOBOO 

II 0000900 

00001 000 

00001100 

00001210 

00001300 

00001100 

00001500 

00001600 

00001700 

00001600 

00001900 

00002000 

00002100 

00002200 

00002300 

00002100 

00002500 

00002600 

00002700 

00002800 

00002000 

0  0  0  0  3  0  0  0 

00003103 

00001200 

00003300 

00003100 

00003500 

OOOOiuOO 

00003700 

O0003HOO 

0000390) 

000010 00 

OC00HO0 

00001200 

00001300 

OOOOIUOO 

0QOOU50O 

00005500 

00005600 

00005700 

000051)00 

0000590  0 

OOOOfeOOO 

00006100 

00O0b20O 

60006500 

00006100 

00  0  06500 

00006600 


en 
en 


f-Sg 


6700 
6600 
6900 
7000 
7100 
7200 
7300 
7100 
7500 
7600 
7700 
7800 
7900 
600ft 
8100 
8200 
8100 
8400 
8500 
8600 
8700 

eeoo 

69no 

9000 

9100 

9200 

qlOO 

9100 

9500 

qhOO 

9  7  00 

9800 

9900 

10000 

10100 

10200 

102IO 

1022O 

102i0 

1025O 

10260 

10270 

10260 

I  0  1  0  0 
10100 
105H0 
10600 
10700 
10600 
I0900 

I I  000 
11100 
11200 
1  1  V10 
1  1100 
11500 

1  uoo 

11700 

naoo 

11900 
12000 


c— 

c— 


OVERRIDE  SOIL  MOIST,  IF  NECESSARY 

IDAY«0 

WRlT£(6,flO0) 
«00  rOaiATC  '  ,  '  DO  YOU  KANT  TO  INPUT  NEW  SOIL  MOISTURE  VALL'E'j 
*'  FOR  BEGINNING'/'  OF  CALC.T  (YESsI  ,NOsO)  IF  NO,  THE  CALC.'i 
*'  WILL  START'/'  AT  THE  TIME  SNOW  IS  OFF  &  «!'H  H20  CONTENT', 
*'  AT  F,C. '/) 

READ{5,720)  10 
720    FOSHAT(lJ) 

IFtIO.EO.0jGo  TO  110 
WRITE(6,1?0> 
120     FOR'UT(/iSx» 'INPUT  AS  DAY  »   AND  THEN  THE  DECIHALCO  TO  l.0)« 
*/'   OF  AVAILABLE  WATER  IN  THE  WHOLE  PROFILE'/) 


730 
«L0 
C— — 
C— .- 


17 


18 
16 


2002 
2003 


11 
11 

12 

11 


C  — 

c— 
c— 


REAn(5,730)  I0AY,SKO 
FORmAT(I1,F3.2) 

Continue 
read  normal  tev1p3 

1U»8 

Continue 

KaJfl 

IF  (IDTYPE.EO.O)  READ t IU, 38,ENDi2002) 

*f  (X'iT(L,N),Nsl,2),XNP<L>/L«0»K> 

IF  (IOTYPE.NE.0)  READ(IU,39iEND»2002) 
*C(XNT(L,N),Nai,2),XNP{L),L*J»K) 
F0Ri1AT(lSX,S(2Fl|.l,F«.2)) 
F0RMAT(15X,5(2F4.1,F«,3)) 

J  =  Jt5 

IF  (J.GT.NOT)  GO  TO  2001 
GO  TO  17 
NDT=L-1 
CO'J  T INUE 
IF(IDTYPE.EQ.O)  GO  TO  11 
Oil  11  I  x a t  ,NDT 

XNPfIX)tXNP(IX)/10. 
CONTINUE 
GO  TO  11 
DO  12  Ix»l.NDT 

X'iP(IX)*XnP(IX)*2.54 
CONTINUE 
CONTINUE 

INPUT  ACTUAL  TEMPS  AND  PRECIP 


WRITE(6,80) 
80  FORMAT ('  ','    00  YOU  HAVE  ANY  ACTUAL  TEMPS  »  PRECIP  TO  INPUT1, 
*'  (YES  =  1,  >IO  =  0)'/t  IF  YES  THEN  ENTER  FIRST  l   LAST  DAY  NUMBERS'/) 
READ(5,fl2)  I&dF.IL 
82  F0R'liT(Ili2lD 

IF  (IO.Ei).O)  NEWTEMbO 
IF  (I.), £0,0)  GO  TO  102 
t»E*TEMsl 
WRlTE(6,81) 
61  FORrUTC  ','  ENTER  VALUES 
*'   OSE  1  COLS 
DO  39  L=IF,IL 
IF  ( I JTVPE.EQ.O)  GO  TO  as 
REAO(5,P7)(xNT(L/K),K=l,2),XNP(L) 


THX.TMN, PRECIP  -(EACH', 
WIThOUT  A  DECIHAL  P0INT)<CR> ' /) 


00006700 
00006800 
00006900 
00007000 
00007100 
00007200 
00007100 
000&7UOQ 
00007500 
00007600 
00007700 
00007300 
00007900 
00008000 
00008100 
00008200 
00008100 

ooooaiao 

00008500 
00008600 
00008700 
00008800 
00006900 
00009000 
00009100 
00009200 
OOOOqlOO 
00009100 
00009500 

onooqhoo 

00009700 
00009300 
000 1)  9900 
00010000 
00010100 
00010200 
00010210 
00010220 
00010210 
00010250 
00010260 
00010270 
00010260 
O001010O 
00010100 
00010500 
00010600 
00010700 

. oooi  onoo 

00010900 
00011000 
00011100 
00O1120O 
0001  1100 

oooi iioo 
ooomiio 

0001 IfaOO 
0001 1700 
OOOI  I  000 

0001  uoo 
00012000 


en 
en 


m 


12100 
12200 
12300 

12<J00 
12500 
12600 
12700 
12600 
12900 
13000 
13100 
13200 
13300 
I  3400 
13500 
13600 
13700 
13600 
13900 
14000 
1110  0 
14200 
14300 
t  U  tl  0  0 
14500 
1<J600 
14700 
14600 
1  49oo 
15C00 
15100 
15200 
15300 
15400 
15500 
15600 
15700 
15'SOO 
15900 
16000 
I  6  1  0  0 
I  j200 
16300 
16400 
16500 
16600 
16700 
16600 
16900 
17000 
17100 
17200 
I  7  3  ['  0 
17*00 
I  7?00 
I  7  0  0  0 
17700 
17  0  0  0 
17900 
16000 
18100 


GO  TO  «9 
85   ft£A0C5,»6)CxNTa,K),k»l,2),XNPCL) 
84  FORMAT<2F«,i,Fo,2) 
67  F0BM*T(2Fa.i,F«,l) 
69  CONTINUE 

IF    (IDTVPE.EO.O)    N«130 
IF    (IOTYPE.HE.0)    h»U5 
00    150    J*liNr>T,5 
00    170    M*J,J+4 
CO    TO    N 
135    IXNP(M)eXNP(MMOOO, 

GO    TO    140 
130    IX^P(H}sXSP(M}*JO«,  • 

140    Co    170    U.'  =  l»2 
170    IxNT(H,L'OsxnT(H,LN)*10, 

»i?ITE(IU,160)t(IXNTCLiK),K«liX)»!XNP{|.)#U«J#J»4) 
160    FOR'UT  ( I  faX,  15I«  J 
150    CONTINUE 

LOCK  IU 
3C2     CONTINUE 
IU  =  13 

oo    continue 
o«    continue 

----      INPUT  DATE  SNOH  OFF 


WRITE  C6,fa0) 
60   FORMAT (■  '.'  INPUT  OAy  #  SNOW  IS  OFF  THE  SITE.  AND', 
*'  THE  nUHBER  OF  DaYj  IN  THE'/1  CROWING  SEASON  FoH  THE',' 
*'  SITE'/) 

PEAD(5,8S9)  ISnOw,NOAVS 
959   FoR'UT(2I3] 

61   CALL  mOIST(IDAY,SMO,NDT,IDTYPE,ILEAP,ISNOW,NDAYS) 


--"  CALC,  GDH  FOR  NEW  TEMPS, 

CALL  r.OH  (IF,  I  L,TL/ TO,  TU,  CONST,  NOSM,NOATM,  I  Y8,ILEAP#NDT»JD0P, 
•IOTYPE) 
00      CONTINUE 
00      CONTINUE 

00  CONTINUE 
20      CONTINUE 

1  CONTINUE 

---«     PREDICT  pHENOLOGICAL  DEVELOPMENT  AND  HTS  OF  SPECIE  SELECTED 


IF  <UROP.Eq.0),0R,(IDOP.EQ.2)> 

*CaLl  PHEN0(I5p> 
LOCK  fl 
LOCK  13 
GO  TO  20 
»999     CONTINUE 

STOP         / 
CMC 
SUBROUTINE  HOlST(IDAY,SHO,IMDT,IDTYPE,IL£AP,ISNOW,NDAYS) 

:-. —      SU6R0UTL  TO  CALC.  SOIL  HOIST  BUDGETS  (EOIV.  DEPTHS) 
:---•      VIA  ET  AND  PCPN  (AFTER  HANKS  1971) 

COMMON  /SLKl/  XNP(475),XNT(475,2),ETT(475),EA(<I75) 
COMMON  /SLK3/  Afc(«75),TRAT(«) 


00012100 
00012200 
00012300 

00012400 
00012500 
00012600 
00012/00 
00012600 
00012900 
00013000 
00013100 
00013200 
00013300 
00013400 
00013500 
0001  JfaOO 
00013700 
00013800 
00013900 
00014000 
00014100 
00014200 
OO0I4J0O 
00014400 
00014500 
00014600 
00014700 
00014800 
00014900 
00O15000 
00015100 
00015200 
00015300 
00015400 
00015500 
00015600 
00015700 
00015300 
00015900 
00016000 
00016100 
00016200 
00016300 
00016400 
00016500 
00016600 
00016700 
00016800 
00016900 
00017000 
0  0017100 
00017200 
00017300 
0  0  017  4  0  0 
00017500 
00017600 
0001 7700 
00017800 
00017900 
00013000 
000)3100 


en 


mi 


18200 
16300 

18100 
16500 
16600 
16700 
166O0 

ie<5oo 

19000 

19100 
19200 
19250 
19300 
19D0O 
195O0 
19600 
19700 
19800 
19900 
20000 
20100 
20200 
20300 
201OO 
205O0 
20600 
20700 
2C3O0 
20900 
21000 
21100 
21200 
21300 
21100 
21500 
21600 
21700 
21000 
21900 
22000 
22050 
22100 
22200 
22300 
221O0 
22500 
22600 
22700 
22800 
22900 
23000 
23100 
23200 
21300 
23-100 
23500 
2JS00 
21700 
23500 
23*00 
2^000 


DIMENSION  PSH(500),LOH(15) 
C  — •- 

C-.-»  SET  LAST  DAY  NUMBERS  OF.  MONTHS 
C . 

DATA  LOM/30, 6 1,91, 122, 153, 181, 212, 212, 273, 303, 331, 36S, 
*J95,«26,«56/ 
C---- 

C— -m  THRESHOLD  LEVEL  FOR  EvAP,  ALGORITHM 


C  — 


5 

10 
c 

c— . 

c— — 

c 


DATA  FCWP/127,,55./,  PTH/5.1/ 

V.CSET=FC 

XH'JTHal.O 
IF  (ILFAP.EQ.O)  GO  TO  10 
DO  S  Is6,15 
LOHtn^LDHCI)*!- 
CONTINUE 
CONTINUE 

SET  WATER  CON,  FOR  FIRST  DAY  EQUAL  TO  FIELD  CAP, 

BS»=wp 

INPUT  NEW  SOIL  MOIST  IF  NEC, 


IF(ir>AY,EO.0)GO  TO  <i00 
PSM{IDAY)s(5HO*(FC-wP))tWP 

BSMcPSH(IDAT) 
CO  TO  U00 

1RITE{6,«10) 
110   ,   F0R1AT(/,5X, 'UPDATE  DATE  IS  POST  DATE  or  PEAK  PftODN,  THUS' 

*/•   PtfOD  EST,  BASED  UPON  TRANS  RATIO  OF  THE  UPDATE  OATCV) 
100     CONTINUE 
C-— 

C— •«  BEGIN  DLY  BUOCETTING  OF  SOIL  HATER 
C 

00  97  HN=1,4 
97  °TRiT(HN)=0.0 

55  =  0, 

56  =  0. 

00  60  J= 1 , INDT 

IF  (J.E0.182)  8SH.FC/2. 

C— —  DET  TIH£  SINCE  LAST  PCPN  GT,  PTH  FOR  EyAP.  ESTIMATION 

THE  =  10. 

11  =  10 

If    (J.EO.l)  Co  TO  63 

IF  ((J.LE.lO).AND.CJ.GT.l))  II«J-1 
DO  61  Ll=l,II 
IL'LL-I 

IF  (LL.EQ.l)  PCPNnXNP(J-IL) 

IF  (LL.GT.l)  PCPNciNP(J«Il.).XNP(J«LL) 
IF(PCpki.GT,PTH)TI«E  =  LL 
IF(p:Pn.GT.PTh)Co  TO  63 
61  CONTINUE 
63  CONTINUE 

c— - 

C---«  DETERMINE  MONTH 
C-  — - 

DO  110  Ilat.ll 

IFfJ.Gf  .i.^H(H).ANO.J.LT.LDM(II  +  l))GO  TO  120 


00018200 
00018JOO 
00016100 
00018SOO 
00010600 
00018700 
00018800 
00018900 
00019000 
00019100 
00019200 
00019250 
00019300 
00019100 
00019500 
00019600 
00019700 
00019SCO 
00019900 
00020000 
00020100 
00020200 
00020300 
00020100 
00020500 
00020600 
00020700 
00020600 
00020900 
00021000 
00021  100 
000212Q0 
00021300 
00021100 
00021500 
00021o00 
00021700 
00021e00 
00021900 
00022000 
00022050 
00022100 
00022200 
00022300 
00122100 
00022500 
00022600 
00022700 
00022*00 
00022 ')  00 
OOf'23000 
00023100 
00123200 
00023200 
000231:00 
00  0  23500 
00023600 
00023700 
O0023eO0 
00023900 
00021000 


CD 

CO 


I 


24100 

24200 

anioo 

24400 
24500 
24600 
2U800 
2u900 
25000 
25100 
25150 
25160 
25170 
25175 
25200 
25100 
25'400 
25500 
25o00 
25700 
25800 
25900 
2b000 
26050 
2bl00  . 
26200 
26100 
26100 
26500 
26600 
26700  ■ 
26300 
26900 
27000 
27025 
27050 
27100 
27200 
27300 
27400 
27500 
27600 
27700 
27800 
27900 
26000 
28100 
28200 
28300 
28400 
26500 
•   28600 
2  8  7  0  0 
26  80  0 
289(10 
29  0  0  0 
29100 
'  29200 
29300 
29^00 
29500 


no 

120 

c— .. 

c— - 
c— - 

c— 

c— - 


64 


c— -• 
c 


70 


c— — 
c 


CONTINUE 
IMsn 

CALC,  POT  ET  FOR  PERIOD 
CALL  TECJiIMJ 
CALC.  ACTUAL  ET 

CALL  ETSEPCBSM,W,TIME,TA,TP,EP,J,lSNOW,NDAYS) 
IF  (J.6T.18H  GO  TO  64 
£A(J)*0.0 
GO  TO  70 
IF  (J.EQ.182)  £TT<J)»0.0 
AW(J)sW 
S5aSS+TA 
S6=S6tTP 

lF(Sfc.EO.O.)st»». 

^CALCULATE  PRECIPITATION  FOR  PERIOD 

IF  tJ.EO.11  SpsxNP(J) 
IF  (J.GT.l)  SP»XNPtJl-XNP(J-t) 
IF  (SP.LT.O.O)  SPoO.O 

PREDICT  SOIL  MOIST  AND  ADJUST  FOR  PHP  AND  rC. 

P5H(J)*BSM-EA(J)-TA+SP 

IF  lPs«(J).CT.FC)  PsMtJ)sFC 

IF  (PSM(J).LT.NP)  PSM(J)orfP 

BS»SPSH(JJ 

GO  TO  60 

IF  (J.EO.ll  SP=XNP(J) 

IF  (J.Gl.lJ  SP-XnP(J)-XnP<J-1) 
IF  (SP.LT.0,0)  SP»0,0 

AW(J>«d 

BS»sPSH<J) 

E*(J)=0.0 
CONTINUE 
RETURN 
END 

SUBROUTINE  TE(LLO.IM) 

Co»*OM  /BJ.KI/  X,NP(«75),XNT(H7S,2j,ETT(fl7S),EA(«75) 
DIMENSION  P0(t5)i  XNDSC15) 


60 


SET  X  OF  DAY  LIGHT  HOURS  PER  MONTH 


OATA  PD/8.J8, 7. 75, 6. 72, 6. 52, 6. 76, 6. 73, 8. 33, 8. 95, 10. 02, 
M 0.08, 10. 22, 9. 51, 8. 38, 7. 75, 6. 72/ 
C  —  -■ 
C SET 

DATA  xNDS/10, 31,30,31, 31. 28, 31, 30, 31, 30, 31, 31, 30, 31, 30/ 


THE  NUMBER  OF  OAyS  PER  MONTH 


C 

C 

C-  — . 


15 


CALC.  AVERAGE  TEMP. 

IF  (IDTYPE.EO.O)  GO  TO  IS 
TAM(tXHT(LLO,n+XNT(LLD,2n/2.)«(">./5.))+32. 

CO  TO  20 
TA=(XMT(LLD,l)+XNT(LLO,2>)/2. 


00024100 
00024200 
00024300 
00024400 
00024500 
00024600 
00021600 
00024900 
00025000 
00025100 
00025150 
00O25160 
00025170 
00025175 
00025200 
000253^0 
00025400 
00025500 
00025600 
00025700 
0002530O 
00025900 
00026000 
00026050 
■  00026100 
00026200 
00026300 
00026400 
00026500 
00026600 
00026700 
00026800 
00026900 
00027003 
00027025 
00C27050 
00027100 
00027200 
00027300 
00027400 
00027500 
0O027b00 
00027700 
00027800 
00027900 
00028000 
00028100 
00023200 
00O28300 
0002B400 
00028500 
00026600 
00028700 
0  0  0  2  8 1!  0  0 
00  0  26900 
00029  0  00 
00029100 
00029200 
00029300 
00029400 
00029500 


en 


29600  20   Er{({TA)*PO(lMn/(l00*XNOS(lM)>) 

27700  10  CONTINUE 

29600  ETTCLLD)=E*25.1 

29900  RETURN 

300PO  END 

30100  SUBROUTINE  ETSEP(BSM,w,  TIME.  TA,  TP,£P,UD#  I3N0H,NDAY5) 

30200  C  —  - 

30J00  C-  —  -CALC.  EA  AND  T*  SEPARATELY  FROM  HANKS  (1971) 

10100  C--  — 

10500  COMMON  /BLK1/  XNP(U75) ,  XNH175,  2)  ,ETT  (175)  ,EA(«75) 

30600  COMMON  /BLKl/  AW(175) , TRAT (1) 

30700  DIMENSION  C0(5> 

30800  C-«»» 

10900  C  — -    SET  COEFFS.  OF  EQN.  OF  EVAP,  FACTOR  VS  TENTHS  OF  GROW,  SEASON 

31000  C-  —  - 

31100  DATA  CO/101. ,-26, 82, -.652, .615, ".0399/ 

11200  ,.»■    DATA  FC.WP/12^.»55./ 

31300  C-— 

Jl<lOQ  C  — --       C*LC.  FRATION  OF  GROW.  SEASON  COHPLETED  AS  OF  TODAY 

31500  C----  ' 

31600  X  =  UO 

11700  YalSNOW 

31000  ZaNOAYS 

31000  F"AC=  CX-Y5/Z 

32000  IF(FRaC.Gt.1.)FBAC«1, 

32100  IFCFSAC,LT,0,)FRAC»0. 

32200  FRAC=FRAC«10. 

32300  C«" 

32100  C--«.    CALC.  WHAT  %   OF  TOTAL  ETP  COULD  GO  TO  POT,  EVAP, 

32500  C---« 

32oOO  FEP  =  CO(l)*-CO(2)*FRACtCO(3)*FfiAC**2,  +  CO(13*FRAC**3.*CO(5) 

'  32700  »»FRaC**«. 

32600  IF<FEP.GT,100.)FEPalOO. 

32900  IF(FEP.LT,0.)FEP=0, 

31000  EP=t(FEP*.01)*ETT(LLD)) 

31100  C---- 

31200  C  —  —   CALC  ACTUAL  EVAP, 

11300     C 

11400  EA(LLD)=EP*CSQRTU/TIME)1 

13500  C-— 

33b00  C.---    CALC.  POTENTIAL  TRANS, 

33700  C— — 

33600  TP=ETT(LLD)-EP 

33900  BB=B5M 

30000  IF  (BH.GT.FC)  BB=FC 

31100  IF  (BB.LT.WP)  BBerfP 

31200  C  —  -» 

31300  C----H  IS  THE  PERCENT  OF  AVAILABLE  WATER 

31100  C-»« 

31500  W=l.-((FC-flB)/(FC-V«P>> 

3u600  C---- 

31700     C CALCULATE  ACTUAL  TRANSPIRTATION 

3uBO0  C---- 

31900  TA=TP»w*2. 

15000  '  IF(TA.CT.TP)TA=TP 

35100  RETURN 

35200  END 

35300  SUBROUTINE  PHENO (NSPECE) 

35100  COMHON  /8L.K2/  |<G  ( 1  75 ) ,  SPK  (2, 1  0 ) ,  nDCDH(175) 

35500  DIMENSION  ACGDM(6,8),HGriTS(6,8),BI0MA5(6»fl)# 

1  r.  h  n  .i  •  I  1  •  v  M  'O 


00029600 
00029700 
00029B00 
00029900 
00030000 
00030100 
00030200 
00030300 
00030100 
OOQJ0500 
00030  600 
00030700 
O0C3OEO0 
00030900 
00031000 
00031100 
00031200 
0  0  011300 
00031100 
00011500 
00031600 
00011700 
00031800 
00011900 
00032000 
00032100 
00032200 
00032300 
00032100 
00032500 
00032600 
00032700 
00032600 
00032900 
00033000 
00033100 
00033200 
00033300 
00031100 
00013500 
00033600 
00033700 
00033600 
00031900 
00031000 
00031100 
00031200 
00031300 
■  OOO3HO0 
00031500 
0  0  0  3 1  b  0  0 
00031700 
00031300 
00031900 
00035000 
00015100 
00035  200 
00035100 
00035100 
OOC3S500 
0003  5600 


en 

o 


— I 


35650 

35700 

35750 

35760 

35300 

15350 

15900 

35950 

16000 

36050 

36100 

36150 

lb200 

36250 

16300 

36150 

16100 

36150 

37100 

37200 

17100 

17400 

17500 

to 

37600 

37  7  00- 

37300 

20 

37900 

25 

38000 

38  100 

36200 

3  310  0 

27 

38100 

33500 

38525 

35 

15510 

33550 

ie600 

10 

18700 

38800 

50 

38900 

55 

39000 

39050 

39100 

56 

39100 

39100 

39500 

39525 

57 

39530 

39550 

3  9  6  0  0 

60 

39700 

39705 

70 

39710 

75 

39715 

39718 

3  9  7  2  0 

76 

39710 

39715 

19710 

39712 

77 

DATA  ACGDH/O. 0,0, 0,0. 0,0. 0,0, 0,0. 0,0. 0,0. 0,2880,0, 

♦  0.0,0. 0,2773.  0,3228.0,  

♦  5030.  0,51 52.0,  3925. 0,7.3«3. 0,5368. 0,6337, 0,6125.0, 

*6109, 0,7006.  0,11 62 1.0, 7500. 0,0. 0,0. 0,0. 0,0.  0,0.0, 
*0. 0,10111. 0,9188. 0,9571. 0,9927. 0,21 05 3. 0,13120.0, 
«t 10 31. 0,1  1622. 0,1 28 02, 0,1 1716. 0,23988, 0,25 376,0, 

*  12675.0, 12675.  0,169  11,  0,11102. 0,26080. 0,27700.0/ 
DATA    HGIlTS/O.  0,2 100, 0,1  100.  0,2600,  0,0  ,0,0.  0,0.0, 

♦  lilOO, 0,2500. 0<59  00. 0,0. 0,0. 0,0, 0,7700, 0,3700,0, 
*9700. 0,0. 0,0. 0,0. 0,0. 0,5000. 0,1, 0,0, 0,0. 0,0. 0,0.0, 
♦6700.0, 0.0, 0.0, 0.0, 0.0, 0,0, 9200. 0,0. 0,0. 0,0. 0,0.0, 
*0, 0,1 1300, 0,0. 0,0, 0,0. 0,0. 0,0. 0,0. 0,0. 0,0, 0,0.0/ 

DATA  QI OH AS/2100.0, 1800.0, 1300. 0,1800, 0,0. 0,0.0, 

•  3900. 0,3100. 0,2700. 0,6600. 0,0. 0,0, 0,5100, 0,1600.0, 

•  1300. 0,11 100. 0,0. 0,0. 0,6600. 0,5600. 0,6200.0, 

*  13200. 0,0. 0,0. 0,7700. 0,6500. 0,8500. 0,0. 0,0, 0,0,0, 

•  8600. 0, 7200, O/lllOOjO, 0,0, 0.0, 0,0, 9100. 0,7900,0, 

♦  11700. 0,0. 0,0. 0,0. 0,10100, 0,6500. 0,0. 0,0. 0,0. 0,0,0/ 

00  20  1=123,301 
00  10  J=l ,5 
I»ATtsJ*J-t 
IF  (NOGDH(IDATE).LT.l)  GO  TO  20 
CONTINUE" 
IBEG=I 
GO  In  25 
CONTINUE 
HSTAGEsl 

00  50  I=IBE6,U57 
NETr,nHs'JC(l)«NC<IBE6> 
IF  (N5TAGE.CT.fi)  GO  TO  55 
IF  UCGOH(nSP£CE,nSTAGE).E<J.0.0)  CO  TO  35 
IF  (NETGOH.GE.ACGOH(NSPECE,NSTAGE>)  GO  TO  40 

GO  TO  50 

NSUGE-JSTAGE  +  1 

IF    (HSTAGE.GE.fi)    Go    TO   55 

GO    TO    27 
I0AY(1,NSTAGE)bI 
NST*CE=hsTaCEM 
CO;iT!*UE 
NSTAGE=l 

DO    70    I=IEEG,«57 
IF    (HSTACE ,GT,8)    GO    TO    75 

IF    (HG*ts<nSPECE.hstaGE),EQ,0,0)    Go   TO  57 
NETHTs  .■|G(I)-NG(IBEG) 
IF    (NETnT.GE.HGHT3(NSPECE,NSTAGE))    GO   TO   60 

GO    TO    70 

NSTAGEsNSTAGEtl 

IF  (N3TA&E.GE.8)  GO  TO  75 

Go  TO  5b 

I0A¥(2,,1STAGE)  =  I  •  ' 

N5TAGE  =  'iSTAGE*l 
COMTIhuE 
NSTAGE=1 

DO  90  I=!BEG,«ST 
IF  (NST»3E ,GT ,8)  GO  TO  100 

If  (alOnAS('«3PECE,NSTAGE).EO.0.0)  GO  TO  77 
NETnIO=t.G£I)-MG(IB£GJ. 
IF  r;ETliI0.GE.BI0HAS(NSPECE,M3TAGE))  GO  TO  60 

GO  TO  90 

NSTAGE^'ISTAGE  +  l 


00035650 

00035700 

00035750 

00035780 

U0035B00 

00035650 

00015900 

00035950 

00036000 

00*036050 

00036100 

00016150 

00016200 

00036250 

00036100 

00016150 

00036100 

0003bl50 

00017100 

00037200 

00017100 

00017100 

00037500 

00037600 

00037700 

00037800 

00037900 

00038000 

00038100 

O0O36200 

00018100 

00013100 

00038500 

00018525 

00038530 

00033550 

00013600 

000  IS  700 

00036800 

00013900 

00019000 

00019050 

000191 00 

00019100 

00039100 

00019500 

00019525 

00019510 

00039550 

0  0  019  6  0  0 

0  0019  7  0  0 

000197C5 

00019710 

00019715 

00019718 

00019720 

00019710 

00039735 

00019710 

00019  712 

00019711 


<T> 


/. 


397U0 

39715 

39750 

397SS 

39760 

39765 

39770 

39775 

39780 

39765 

39790 

39795 

39600 

39605 

39810 

39615 

39620 

39325 

39030 

39335 

39810 

398i5 

39650 

3990O 

10000 

10100  ■ 

10200 

10300 

10100 

10500 

106OO 

10700 

10800 

10900 

11000 

11100 

nauo 

11300 
11100 
11500 
11600 
11700 
11600 
11900 
12  0  0  0 
12100 
12200 
12310 
12100 
12500 
12o00 
12  7  0  0 
12800 
12900 
1  i  0  0  0 
13100 
13200 
13300 
13100 
13500 
13600 


80 

90 
100 

no 


120 
130 


c— - 
c— - 

c— - 

C-— 


GO   TO   76 

IDAY<3,NSTAG£)«I 

nsuge.snsta&e  +  1 

continue 

FORHAtC"  ",8X,»DaY  «of  phenol,  stage  oevel.v, 

*"  ot  >    2   ,  3    «    5    6  7  e»/JX,IJ,7C2X.lJ>) 

HPnC(7,120)tIOAY(2iJ),Jsl,8) 

FORMAT!"  «,8X,»0aY  *  OF  HEIGHT  LEVELS  (CM)"/, 

*"     5    1°    15    20    2S    30  35  10"/3X,  I3,7(2X,  I  J) ) 

KRIUC7,130)UDAY<3,J),j=l,e)  '* 

FORMATC"  ",8*,"DaY  *  OF  BIOHASS  OEVEL.  CGHs)b/, 

*RETURN    10    '5    2°    "    3°  3S  ao"'3X,I3,7C2X,n)> 
END 

SUBROUTINE  GDH( IF, IL,TL, TO, TU, CONST, NOSH, NOATM.IYB, 
*lLEAP,NDT,IOop.It>TYPE) 

CQ"".ON  /BLKly  XNP(175),XNT(175,2),ETT(«7S),EA(47S) 
COMMON  /iiLK2/  NG(175),5PKt2,10),NOGOHCl75) 

CO'HOH  /SLK3/  A-.C175),TRATU)  i 

DIMENSION  NDS(15),XHN(1S)  ' 

DIMENSION  TH(2l),RH<21) 

OATA  "Cth/0,25/,YInT/0.5/ 

DATA  NOS/ 30, 3 1,30, 3 1,3 1, 28, 3 1,30, 3 1,30, 31, 3 1,30, J  1,30/ 

DATA  x.i'i/'SEP','OCT',  I  NOV'.  'DEC,  '  JAN  I ,  '  FEB ' ,  «  HAR  • ,  •  APRI  . 
••may1, 'JUN», "Jul1, 'AUG1, 'SEP1, 'OCT', 'NOV'/ 

DATA  NP/12/,Nl/2/,OTN/|7.0/,DTX/22,/,FX/.25/ 


SOt 


23 

21 


LOOP  FOR  OAyS 

IY»=IYR 

Sf'GDH  =  0,0 
IF  (ILEAP.EO.l)  NDS(6)*29 
N0H=NOS(11 
I  YE  =  I Y0+ 1 

JMal 
DO  20  1=1, NOT 

IF  (IDTYPE.NE.O)  Go  TO  801 
XNT(I,l)=(XNT(I,l)-32,)*(S./«,) 
XNT(I,2)=(XNT(I,2)-32.)*t5,/9.) 
TAVE=CXHT(I,1)+XNT(1,2))*0.S 
IF  (LEO.])  GO  TO  2J 
Tl=XNT(I-l,2) 
GO  TO  2a 
Tl=XNTfI,2) 
T2=XNT(I,1) 
T3=XNTtI,2) 

C*LC.  HRLY  TEHPS, 


IP 
J  =  NP 


70 


K»l 

Khrll 

XXN'P 

!l  =  l 

DO  50  J=1I,NH 

Z=J 

IF(K 

TriCJ 

RH(J 


,ri].2)ZnJ-NP 
Z-l  .)/(XXNP-l.) 
)=T2-C(T2-T1)«FC> 
)  =  J.O 


"  1 


000397«a 

00039715 

00039750 

00039755 

00039760 

00039765 

00039770 

0003977S 

00039780 

00039785 

000'39790 

00039795 

00O39300 

00039605 

00039810 

00039615 

00039820 

00039325 

00039330 

00039.635 

000J981O 

00039815 

00039850 

00039900 

00010000 

00010100 

OOOU0200 

00010300 

00010100 

00010500 

00010600 

00090700 

00010800 

OOOIO9OO 

00011000 

00011100 

00011200 

00011300. 

00011100 

0001ISOO 

00011600 
0  0  0117  0  0 
000<4 1600 
00011900 
00012000 
00O121O0 
00012200 
00012300 
00012400 
0001*2500 
00092600 
00OU27O0 
0  0  012300 
0  0  012900 
00043000 
00013100 
00013200 
00013300 
00013UOO 
00013500 
00013600 


«J700 

43690 

43900 

44000 

tJ  100 

44200 

4al00 

44400 

4u500 

44600 

V4  700 

46800 

4«900 

0^0  0  0 

45100 

45200 

45300 

45400 

45500 

456O0 

45700 

45800 

45900 

46000 

46tU0 

4  6  2  0  0 

46300 

46400 

46500 

46600 

46700 

46800 

46900 

47000 

47100 

47200 

4?J00 

47400 

47500 

47600 

47700 

47300 

47900 

48000 

46100 

43150 

48200 

48100 

MS400 

46500 

46o)0 

4  8  7  0  0 

48800 

4  8900 

49000 

49100 

49200 

49300 

49-00 

q  9  5  0  0 

H<iDPO 


c». 
c— 

c— 


IF  (THCJ).LT.THI)  GO  TO  SO 

DT»TM(J)-n 

IF  (OT.GE.DTX)  RHCJ)aFX 

IF  ((nT.GE.DTX).OR.(OT.LE,t>TN))  CO  TO  50 

RH(J>=J,55t<(Cl.0-FX)-/CDTN-D.TX)>*DT) 

50  CONTINUE 

If  tK.'lE-.  1)60  TO  60   . 
II=NPtl 

NH  =  24 

XXNP=24-NP 
K  =  2 
T1=T3 
CO  TO  70 
60  CONTINUE 
20S  CONTINUE 

CALC.  COM  FOR  OAY 

L»I 

G0=0.0 
DO  21  K=l,2tt 
FSMsl.O 

fact*!. a 

ArTO-Tl. 

IF  ((JH.OE.1).*MO.(JH.LE.5>)  CO  TO  JO 
IF  (MOSM.EO,0)  Go  TO  10 
IF  tA*'(I).ST.«*CTH)  GO  TO  30 
FSMsYINT*(((ltO-VINT)/WCTH)*AH(Xl) 
10   IF  (NOAT^.EQ.O)  GO  TO  IS 
IF  (RH(X).GT.TL)  GO  TO  10 
A  =  A*(  (FSH«RH(K))/CONST) 

Gn  TO  30 
15  AsA.FSH 
JO  CONTINUE 

TnaTh(K) 

IF  (TN.LT.TL)  TNML 

IF  (IM.GT.TU)  TNaTO 

IF  (TN.GE.TO)  GO  TO  400 

DGD  =  U/2.>Ml'CosU.l«159+(J,i«159«(<TN"Tt.>mo-Tl>)m 

(100  DGoI(A)MlKO5((J.l«15')/2.)  +  ((J.t'tl5V2.)*((TN-T0>/(TU-TO5))j) 

"JIOl  C0  =  G0tOG0 
21       CONTINUE 

SMGP:-(=SKGnH+GO 
NG(I)=S"GOH 

H06DHtI)aGD 
IF  (I.LE.NDhj  GO  TO  6 

JKsJM+l 
ND^  =  'I0'1  +  NDS(JM> 
6  IF  (I.E'J.123)  IYR=I YE  _  . 

IDC:1-I.D'J  +  .',D3U1) 

U  (I.NE.l)  GO  TO  498 
IF  <(IDOP.EO.l).0R.U0OP.EQ.2>>  WRITE  (7#  499) 
q99   PORHATf  >,18X,'TE"PERATDSE'>6Xf  "HEAT  UNI  TS'  ,  6X.  '  EVAPQRA  •  , 

.iTION  UNlTS'/lX.tOAY     DATE      TMX    T«N   TAVE    GOH      SHI, 
»'COH   PKECIP       ETT      EA      A.H20'/) 
<|98   IF  ((lOoP.E0.1).OR.(inoP.EO,2))  WRITE<7,500)  .,-,.. 

.I,IDC,XMN(jM,,IYR,(XNT(I,L),Lal,2),TAVE,GD,SMG0H,XNP(n.ETT(n» 

500  FORMAT l"    ',I4,l3,A3,I4,2F6.1,F6,2(F7.2,F9,l.tX,F5.J,JF9.«) 

20  COHTI'.UE 


oooaJ70fl 

00043800 
00041900 

00044000 

00044  ion 

00044200 

00044300 

00044400 

00044500 

00(144600 

00044700 

00044BO0 

00044900 

00045000 

00045100 

00045200 

00045100 

00045400 

00045500 

00045600 

00045700 

00045800 

00045900 

00046  0  00 

00046100 

00046200 

00046300 

00046400 

00046500 

00046600 

00046700 

00046600 

00046900 

00047COO 

00047100 

OC047200 

00047300 

00047400 

00047500 

0004H00 

00047700 

00047600 

00047900 

000480  0  0 

00048100 

00048150 

00048200 

00048300 

00048400 

00046500 

00048600 

0CO4e7O0 

00048800 

00048900 

00049000 

0  00  4  9  10  0 

00049200 

00049300 

00049UOO 

00  0  49500 

00049600 


CO 


64 


o  o  o 

O  <3  o 
MOO- 

c  o-  a 

<T  =»  ■=* 
o  o  o 
o  o  o 


a: 

»- 
ui 
oc  o 

z 


o  o  -o 
o  o  o 
t*-  ff>  w 
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